STOCHASTIC BIFURCATION MODELS 



Richard F. Bass and Krzysztof Burdzy 

Abstract. We study an ordinary differential equation controlled by a stochastic process. 
We present results on existence and uniqueness of solutions, on associated local times (Trot- 
ter and Ray-Knight theorems), and on time and direction of bifurcation. A relationship 
with Lipschitz approximations to Brownian paths is also discussed. 
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1. Introduction. 

Let B t be a continuous function of t, let to, xo, Pi, fa £ R, and consider the ordinary 
differential equation 



A ifX t <B t , 
fa ifX t >B t , 



t e M, 



A" (t ) = zo- 



(1.1) 



Among the results we prove are the following: 

(1) Although in general there will not be a unique solution to (1.1), there will be a 
unique Lipschitz solution to (1.1) if St is a typical Brownian motion path. 

(2) Let B t be a Brownian motion with B = and let denote the solution to 
(1.1) when t = and X(t ) = x . The map y — > X\ is a one-to-one map of H. onto HL 
The smoothness of this map is controlled by the local time at of X\ — B t . If we call this 
local time L\ and fa-,fa satisfy suitable assumptions, then L\ is jointly continuous in y 
and t and {L^, y > 0} and {L~^v , y > 0} are strong Markov processes. We show that this 
implies that for a fixed t > 0, the function y — > A" t y is of class C 1+7 with 7 < 1/2, but it 



(3) As we shall see below, (1.1) is an example of a bifurcation model; if B t is 
a Brownian motion, (3\ < and /?2 > 0, each of the events {lim^oo X t = +00} and 
{lim^oo Xt = —00} has positive probability. The bifurcation time is defined by = sup{t : 
X t = B t }. We calculate both the probability of {lim^oo X t = +00} and the expectation 
of the bifurcation time using excursion theory. 

(4) The equation (1.1) sheds light on the best Lipschitz approximation to Brownian 
paths. In particular we obtain an estimate on the lower bound on the best constant in the 
Komlos-Major-Tusnady result concerning strong approximations of Brownian motion by 
random walks. 

Equation (1.1) is similar to an equation that arose in the course of an economic 
study and its accompanying probabilistic model in Burdzy, Frankel, and Pauzner (1997, 
1998). These papers introduce and study an economics model whose technical side is based 
on the following equation: 



where B t is a Brownian motion starting from B = bo, [3 > is a fixed constant, and / is a 
non-increasing Lipschitz function. The case when xq = /(&o) is of special interest. Results 



is not C 3 / 2 . 




t > 0, 




(1.2) 
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on the time and direction of the stochastic bifurcation were crucial elements of these two 
papers. 

We also consider the following equation, more general than (1.1). 



dX t __ 


f Pi 


x t 


-B t \ ai 


if X t 


< B t , 


dt 


\p2 


x t 


-B t \ a * 


XX t 


>B t , 



t e 



X(t ) = Xq. 



(1.3) 



(If cti = «2 = 0, then (1.3) reduces to (1.1).) Equation (1.3) was inspired by the following 
model. Consider a pendulum with rigid arm which is turned upside down (see Fig. 1.1). 



w 




A 

Figure 1.1. 



Let X t denote the distance of the weight W from its unstable rest position at the 
top of the vertical arm. When X t = x and x is small, the weight is about c\x 2 units below 
its rest position and, therefore c^x 1 units of potential energy must have been converted to 
kinetic energy, given by c^(dX / dt) . Hence, we have the approximate relationship dX/dt = 
c^Xt, assuming infinitesimally small velocity at the rest position. Note that if the initial 
velocity at the rest position is close to zero, then the time it takes the pendulum to move 
any fixed non-zero distance from the rest position is very large. We now add stochastic 
oscillations to our pendulum model. We suppose that the base A of the pendulum vibrates 
according to a Brownian motion B t . Then the position X t of the weight W relative to 
A is X t — B t and we have dX/dt = c^{X t — B t ), which is (1.3) with a\ = ai = 1 and 
-P 1 = p2 = c 4 . 
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The solutions to (1.1) exhibit fast switching between two kinds of excursions. See 
Karatzas and Shreve (1988, Sect. 6.5) for a closely related model. Mandelbaum, Shepp, 
and Vanderbei (1990) also consider a model with fast switching between two kinds of 
excursions, but we were not able to find a direct connection with our own model. 

The rest of the paper consists of five sections. Section 2 contains results on existence 
and uniqueness of solutions to (1.1), (1.3), and related equations. The process B t will 
generally be a Brownian motion, but Theorems 2.3 and 2.4 also apply to some fractional 
Brownian motions (see Examples 2.3 and 2.4). 

Let Xf denote the solution to (1.1) with Xq = y. For a fixed t > 0, the function 
y — > X\ is a transformation of R onto itself. How smooth is this map? How many 
derivatives does the function y — > Xf have and are they continuous? To answer these 
questions, one is led to study the local time of X\ — B t . Section 3 is devoted to a number 
of results about local times related to (1.1), including analogues of the Trotter and Ray- 
Knight theorems. See Knight (1981), Leuridan (1998), Norris, Rogers and Williams (1987), 
Revuz and Yor (1991) and Yor (1997) for old and new variants of the Ray-Knight theorem. 
Our local times are defined as local times at points, but they may also be viewed as local 
times of Brownian motion on a random curve — see (5.15) in Follmer, Protter, and Shiryaev 
(1995) for a result on local times on non-random curves. 

Section 4 gives explicit formulae for the probability of upward bifurcation for the 
equation (1.3) and the expected bifurcation time for (1.1), with some indication how to 
proceed in the more general case (1.3). This extends results from Burdzy, Frankel and 
Pauzner (1998). Section 5 takes a look at the solutions to (1.1) as Lipschitz approximations 
to the Brownian path. As a consequence we obtain some lower bounds related to the 
Komlos-Major-Tusnady construction; see Theorem 5.7. Finally, Section 6 is a list of open 
problems. 

In Sections 3-5, we consider Brownian motion defined on the whole real line R, 
i.e., the process {B t , — oo < t < oo}, where {B t ,t G (0, oo)} and {B_ t ,t G (0, oo)} are 
independent Brownian motions starting from with variance EBf = EB 2 _ t = a 2 t. Unless 
stated otherwise, we will assume that all Brownian motions (including those with drift 
and/or reflection) have infinitesimal variance a 2 , and that all constants are strictly positive 
and finite. 

Section 3 of the paper was inspired by unpublished heuristic calculations involving 
local times which were a part of an earlier project of David Frankel, Ady Pauzner, and 
the second author. We would like to thank the many colleagues who kindly gave us advice 
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on various aspects of the model: Robert Adler, Ludwig Arnold, Jean Bertoin, Miklos 
Csorgo, Burgess Davis, Laurent Decreusefond, David Frankel, Mike Harrison, Haya Kaspi, 
Frank Knight, Jim Kuelbs, Avi Mandelbaum, Ady Pauzner, Jim Pitman, Philip Protter, 
Emmanuel Rio, Ruth Williams, Marc Yor, and Ofer Zeitouni. 



2. Existence and uniqueness of solutions. In this section we present several theorems 
on the existence and uniqueness of solutions to differential equations similar to (1.1). There 
is considerable overlap among the theorems, but each contains cases not covered by the 
other. We first present our main results. They are followed by some remarks and examples. 
The proofs are relegated to the end of the section. 
We start with the equation 



dX 




\Xt 


-B t \ 


~dt ~ ' 


I ft 


\x t 


-B t \ 



tem, x(t ) = x . (2.1) 



where B t is a Brownian motion, oi, 02 > — 1, and Pi, @2 £ 

First note that the function X t = B t is a solution to (2.1) with to = and xq = 0, 
because neither of the conditions on the right hand side of (2.1) is ever satisfied. We 
would like to disregard such a solution for two reasons. First, the economics model behind 
(1.2) required that the solutions to (1.2) be Lipschitz. Second, the example X t = B t is 
rather artificial. For oi, 02 > it is natural to require that X t is a Lipschitz function. We 
generalize this to all 01,02 > —1 by writing an integrated version of (2.1), namely, 



X t = x + [ [Pi\X s -B s \^l {Xs _ Bsm +p 2 \X s -B s \ a2 l {Xs _ Bs>0} ]ds. 

J to 



(2.2) 



It is easy to see that solutions to (2.2) satisfy (2.1), but the example X t = B t shows that 
the opposite statement is not true. 



Theorem 2.1. For fixed to, xq, ft G R, a 2 > 0, and 0:1,0:2 > —1, there exist a 
Brownian motion B t and a process X t which satisfy (2.2) with the initial condition as in 

(2.1) . The solution X t is unique in law. We may construct X t in such a way that {X t , B t ) 
is a strong Markov process relative to the appropriate filtration. If we assume in addition 
that 01,0:2 > 0, then for a given Brownian motion B t there exists a unique solution to 

(2.2) , a.s. 
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Our next theorem is a result on existence. We will state the result for the following 
generalization of the equation (1.1), 



Theorem 2.2. Assume that Fi and F 2 are continuous functions and that \Fi\ and \F 2 \ 
are bounded by (3 < oo. If B t is a continuous process, then (2.3) has a Lipschitz solution, 
a.s. There exists a maximal Lipschitz solution {X^,t > to} to (2.3); it is adapted to the 
filtration T t = cr(B s , s G [to, £])• 

Haya Kaspi pointed out to us that measurability of a solution to (2.3) is the most 
delicate point of Theorem 2.2. 

We will say that L x is a local time for a process B t if it is the occupation time 
density: 



for all h bounded and measurable. Note that if B t is continuous and the local time Lf is 
jointly continuous, then sup x L x < oo, a.s. for each t. 

We will use the traditional Markovian notation P x to denote the distribution of 
{Bt, t > to} conditioned by {B to = x}, even though we do not assume the Markov property 
for B t in Theorems 2.3 and 2.4 below. 

Theorem 2.3. Let to > 0, xo>/3i>/?2 £ Assume that 

(i) the process B t is continuous and has a jointly continuous local time Lf, and 

(ii) if A t is an adapted process with A to = xq whose paths are Lipschitz continuous with 
Lipschitz constant M, then for each x the law of {B t + A t ,t < t < t + s} under ¥ x 
is mutually absolutely continuous with respect to the law of {B t , t < t < t + s} 
under P x + x o ; for every s > 0. 

Then with probability one there exists a random so > and a unique Lipschitz solution 
to (1.1) on [t , t + sq\. 

If in addition we assume that B t is strong Markov then there is a unique Lipschitz 
solution to (1.1) for all t >t . 

Remark 2.4. If W t is a Brownian motion and / is a strictly increasing function such that 
both / and / _1 are Lipschitz continuous, it is easy to check that B t = f(W t ) is a strong 
Markov process that satisfies the other assumptions of Theorem 2.3. 





(2.3) 




a.s., 
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Theorem 2.5. Let t , x G M. Assume that F\ and F 2 are bounded, Lipschitz functions. 
Suppose that both are bounded by M and that both have Lipschitz constant less than or 
equal to M. Let B t be a continuous process such that 

(i) there exist c\ > and 7 G (0, 1) such that whenever s < t, 

¥(B t G dy I F a ) < -g— dy, y G R, (2.4) 

(ii) if At is an adapted process with At = xq whose paths are Lipschitz continuous with 
Lipschitz constant M , then for each x the law of {Bt +At, to < t < to + s} under F x 
is mutually absolutely continuous with respect to the law of {B t , t < t < to + s} 
under p x + x o ; for every s > 0. 

Then with probability one, there exists a unique solution to (2.3) for all t > t . 

We will show in Example 2.10 below that Theorem 2.5 applies to some fractional 
Brownian motions. As in Remark 2.4, some functions of fractional Brownian motions also 
satisfy the hypotheses of Theorem 2.5. 

Let f(x,b) = f3il{ x <b} + (^2^{x>b} and suppose that a\ = ol<i = 0. Then (2.2) may 
be written as 

X t = x + [ f(X t ,B t )ds. (2.5) 

Jt 

The function (x,b) — > f(x,b) is discontinuous. In applications, such as that in Burdzy, 
Frankel, and Pauzner (1997), it may be argued that a model with continuous dX/dt might 
be more realistic. Let us replace / with a continuous approximation, 



f e (x, b) = Pil{ x<b - e } + /3 2 l{x>b+e} + 

and consider the corresponding equation 



P2~P 



2e 



-(x -b + e)+p! 



l{f>-e<x<6+e}7 



Xf = x + [ f s (Xf,B t )ds. (2.6) 

J t 

We will show that the solutions to (2.6) converge to those of (2.5), and thus many results 
about solutions to (2.5) proved later in this article may be applied to give asymptotic 
results for the solutions to (2.6). 

Theorem 2.6. Assume that the equations (2.5) and (2.6) are defined relative to the same 
Brownian motion B t . The equation (2.6) has a unique Lipschitz solution. As e — > 0, the 
functions Xf converge to the unique solution X t of (2.5), a.s. 
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Note that the convergence in Theorem 2.6 is uniform on compact sets as all functions 
Xf are Lipschitz with constant max{|/?i|, \P 2 \}. 

Remark 2.7. For the economics model behind (1.2), one does not necessarily want to 
require the Markov property to hold. The proof of Theorem 2.3 uses the strong Markov 
property to do an induction argument. For Theorem 2.5 we have in mind examples where 
B t is a Gaussian process; see Example 2.10 below. In general, B T+t — B T will not be 
Gaussian when T is a stopping time. 

Example 2.8. We present an elementary example of a continuous deterministic function 
t — > B t for which there are multiple solutions to (1.1). Let (3\ < 0, p 2 > 0, 

((l + p 2 )t forte [0,1], 
B t = <l + /3 2 for t > 1, 
I for t < 0. 

There are uncountably many solutions to (1.1) with this choice of B t and the initial con- 
dition Xq = 1. Here are two of them: 



for t < -l//3 2 , 

l+(3 2 t for£>-l//? 2 ; 



n = 



for t < -l/p 2 , 
_ J 1 + fot forte (-1/02,1] 

1 + (3 2 forte (1,5], 
1 + (3 2 + 5/3i + Pit for t > 5. 



Example 2.9. As we noted earlier in this section, X t = B t is a solution to (1.1) but a 
rather trivial one. In this example, we will show a less trivial and perhaps more interesting 
non-Lipschitz solution to (1.1). Take Pi = p 2 = in (1.1); in other words, consider the 
equation 

^ = if X t ^B t teR, X(t ) = x . 

The function X t = is a solution to this equation and, moreover, it is the only Lipschitz 
solution, by Theorem 2.1. Let Y t be a skew Brownian motion, i.e., a process which may 
be constructed by flipping positive excursions of a standard Brownian motion B t to the 
negative side with probability p\ and negative excursions to the positive side with proba- 
bility p 2: independently of each other. Suppose that pi ^ p 2 so that the process Y t is not a 
standard Brownian motion. Let L t be the local time of Y t at 0. By a result of Harrison and 
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Shepp (1981) (see also Exercise X (2.24) in Revuz and Yor (1991)), for a suitable constant 
ci 7^ 0, the process Y t — c\L t is a standard Brownian motion. If we take B t = Y t — c\L t 
then X t = c\L t is a non-Lipschitz solution to our equation. 



Example 2.10. We provide an example of a process satisfying the assumptions of Theorem 
2.5 that is not strong Markov. Let B t be fractional Brownian motion of index H E (0, 1/2]. 
This means that B t is a mean zero Gaussian process with 

Cov(B s ,B t ) = Cl (s 2H + t 2H -\t-s\ 2H ). 

B t has a stochastic integral representation 

B t = I R{t,u)dZ u , 

J — oo 

where Z u is a standard Brownian motion and 

R(t,u)=c 2 [((t-u) + ) H - 1/2 -(u-) H - 1 / 2 ]; 

see, e.g., Rogers (1997). Conditioning on T s with s > 0, the law of B t given T s is that of 
a Gaussian process with variance 



4E 



(t-u) H -^ 2 dZ u y \F 8 ] =c 2 £ (t-u) 2H - 1 du = c 3 (t-s) 2H . 



Assumption (i) of Theorem 2.5 is immediate from this. 

We now show (ii). We give the argument for the case t = xo = 0, s = 1; the 
extension to the general case is routine. 

If H = 1/2, then B t is standard Brownian motion, and (ii) follows from the Girsanov 
theorem; so we suppose H < 1/2. Let a = H + 1/2. See Decreusefond and Ustiinel (1997) 
for more details of some of the steps in the following argument. Let F(a, 6, c, z) be the 
standard Gauss hypergeometric function and define an operator Kh on functions on [0, 1] 
by 

(K H f)(t) 

= r(g + l/2) I {t " X ^~ 1/2F ( H ~ X / 2 ' l / 2 ~ H ' H + V2, 1 - t/x)f(x)dx. 
Let H H = {K H h : h e L 2 ([0, 1])} and define 

H H = II^hVIIl 2 - 
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For (3 E (0, 1) define 

(I P f)(x) = f ^- ) J Q X f(t)(x-tf- 1 dt 

and 

By Decreusefond and Ustiinel (1997) (Theorem 2.1, Theorem 3.3, and the proof of Theorem 
3.3), we have that Tin is dense in the set of continuous functions on [0, 1] that are null at 
and that Kh is an isomorphism from L 2 ([0, 1]) onto / H+1 / 2 (L 2 ([0, 1])). By Proposition 
2.1 of that paper, is the inverse to I 13 . 

Since K^ 1 is continuous from I H+1 / 2 (L 2 ) into L 2 , then KJj 1 oI H+1 / 2 is continuous 
from I? into itself, and so there exists a constant c 4 such that 

\\K^I H+1 l 2 g\\ L ^ <C4|M| L3 . 

Thus if / e H H , then 

\\K^f\\ L 2 <c4D H+1 ' 2 f\\ L ^ 

or 

||/||w H <c 4 \\D<*f\\ L ,. 

Let A t be a uniformly Lipschitz process as in the statement of Theorem 2.5. By 
Theorem 4.9 of Decreusefond and Ustiinel (1997) and the Novikov condition discussed just 
after that theorem, (ii) will hold if for each T e (0, 1) we have 

Eexp[P(.)||^/2]<oo. 

By the above paragraph, it is enough to show 

Eexp(jf \D a At\ 2 dt/2^j < oo. (2.7) 

To show (2.7), by an approximation argument it suffices to show that for each fixed T > 
there exists c 5 (depending on T) such that if / is a C°° function on [0, oo) with /(0) = 0, 
then 

sup \D a f(x)\ <c 5 ||/'||oo; (2-8) 

0<t<T 

(2.7) will then follow easily from (2.8) and our assumptions on A t . 
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Note that by a change of variables, 

I 1 - a f(x)=c 6 [ f(x-t)t- a dt, 
Jo 

and by the Leibniz formula and the fact that /(0) = 0, 

^ 1_c 7(z) = c 6 £ f'(x - t)t~ a dt = c 6 £ f'(t)(x - t)~ a dt = I l ~ a f'(x). 

Since a = H + 1/2 < 1, then \x — t\~ a is integrable on [0, x]. So, for u = /', 

\D a f(x)\ = ll^uix)] < H^Hoo / \x-t\~ a dt < crllulloo 

for x < T. This gives (2.8), and thus a fractional Brownian motion with parameter 
-H" e (0, 1/2] satisfies the assumptions of Theorem 2.5. 

Example 2.11. The weaker version of Theorem 2.3, i.e., the one without the assumption 
on the Markov character of B t , applies to fractional Brownian motions with parameter 
H E (0,1/2]. Assumption (ii) of Theorem 2.3 is the same as (ii) of Theorem 2.5; we 
have verified that assumption in the previous example. As for assumption (i) of Theorem 
2.3, the joint continuity of the local time for the fractional Brownian motion follows from 
Lemma 8.8.1, Theorem 8.8.2 and the proof of Theorem 8.8.4 in Adler (1981). 

Example 2.12. Fabes and Kenig (1981) gave an example of a process B t satisfying 

dB t = a(B t ,t) dW t , 

where Wt is a standard Brownian motion, a is Holder continuous in the first variable, a is 
bounded above and below by positive constants, and the distribution of Bi does not have 
a density with respect to Lebesgue measure. B t is a space-time strong Markov process. 
Because a is bounded below, it is not hard to see that B t has a jointly continuous local 
time (cf. Revuz and Yor (1991), Ch. 6) and that hypothesis (ii) of Theorem 2.3 holds. 
Thus this process B t is an example where the assumptions of Theorem 2.3 hold, but those 
of Theorem 2.5 do not. 

The rest of the section contains proofs of our main results. The following lemma is 
immediate. 
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Lemma 2.13. 


Let B t 


= B. 


-t and 


x t 




X. 


-t- 


solution to 
















dX 


I -Pi 


\Xt~ 


B t \ ai 


if 


x t 


< 


B t , 




I -ft 


\x t - 


B t \ a * 


if 


x t 


> 


Bt, 



t e R, X(-t ) = x . 



Proof of Theorem 2.1. For simplicity, assume that to = 0. The equation 

Y t = x + [ [p 1 \Y s \ a n {Ysm + p 2 \Y s \ a n {Ys>0} ]ds- I dB s , t>0, 
Jo Jo 

has a weak solution which is unique in law by Theorem 5.15 in Karatzas and Shreve 
(1988). For X t = Y t + B t , the last equation is equivalent to (2.2) for t > 0. This proves 
the first assertion of the theorem. The strong uniqueness in the case ai,a 2 > follows 
from Proposition 5.17 of Karatzas and Shreve (1988). We note that although the function 
y — > y a is not bounded, that proposition clearly applies by using a truncation argument. 
The part of the solution to (2.1) for t < t = can be obtained in a similar way using 
Lemma 2.13. That X t may be constructed so that (X t , B t ) is a strong Markov process 
follows from the weak uniqueness in a standard manner; see Bass (1997), Section 1.5, or 
Stroock and Varadhan (1979), Chapter 6. □ 



Proof of Theorem 2.2. We start by showing that for each uj and for any u\ and z\ there 
exists a maximal solution X^ 1,Zl to the equation 

dX t /dt = F 1 (X t ), teR, X(u 1 )=z 1 . 

First of all, it is well known that there exists at least one solution to the equation since F\ 
is continuous. Since |Fi| is bounded by /?, all solutions are Lipschitz with constant j3 and 
so their supremum X™ 1,Zl is also a Lipschitz function with constant (3. Next note that the 
maximum of any two solutions is also a solution to the equation. This and the Lipschitz 
property of solutions easily imply that there exists a sequence of solutions converging to 
X^ 1,Zl , uniformly on compact intervals. Now a standard argument can be used to show 
that X^ 1,Zl is a solution to the equation. 

The analogous maximal solution to dX t /dt = F 2 (X t ) with the initial condition 
X{ Ul ) = zi will be denoted X" 1 '* 1 . 

We start by proving the existence of a solution to (2.3) for t > to. Consider a small 
5 > 0. We proceed to define a 5-approximate solution Xf to (2.3). First suppose that 
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B to < xq. By the continuity of the paths of B t , for almost every path of B t , there exist 
a unique time t\ G (to, oo] and a function defined for t G (to,t\), such that X^ Q = xo, 
X 5 ti = B tl if t x < oo, and X 5 t = Xl"' x ° for all t G (t , h). We then let X* = X 5 ti + (3(t - h) 
for all t G [ti, ti + 5], if ti < oo. If B to > x we use the same procedure to define Xf for 
t G \pQit\ + 6} except that we use the function xl°' x ° in place of Xl° ,Xi) . If B to = xo, we 
let U = t and Xf = X 5 ti + (3{t - t x ) for t G [h, h+5]. 

We have defined Xf on an interval [to,t\ + 5]. Let xi = X% +5 . Let us replace 
the initial condition in (2.3) by X(t\ + 5) = x\ and define an approximate solution X^ to 
(2.3) on an interval \t\ + 5, t 2 + 5] using the same method as above. By induction, we can 
construct a (possibly infinite) sequence of times {tk} and a continuous function Xj* which 
satisfies (2.3) on every interval (tk + 5, tk+i) and which is linear on every interval [tk, tk + S] , 
for k > 1. Note that the function X^ is defined for all t > to because tk+i > tk + 5 for 
every k. 

By construction, the 5-approximate solution Xf is a Lipschitz function with Lips- 
chitz constant (3. 

For every integer m > 1 consider a 1/m-approximate solution x\^ m . All of these 
functions are Lipschitz with the same constant (3, and they all satisfy X^™ = xq- Let X t 
be defined by 

X t = limsup X t 1//m = lim sup x\^ m . 

The supremum of an arbitrary family of Lipschitz functions with constant (3 is 
a Lipschitz function with the same constant, and the same remark applies to the limit 
of a sequence of such functions. Hence, for every n, the function Y t n = sup m>n x\^ m is 
Lipschitz with constant /?, and the same is true of X t . Note that converge in a monotone 
way to X t , uniformly on compact intervals, because all these functions are Lipschitz with 
the same constant (3. 

We will show that X t is a solution to (2.3). Let 

W{8) = |J {(t, y):y = X + (t- s)(3, te[s,S + 5]}. 

{{s,x):s>t ,B s =x} 

For 5 < 8i, the portion of the graph of Xf which lies outside W(5i) satisfies (2.3), by 
construction. 

The set of t such that B t = X t is closed because both functions B t and X t are 
continuous. Consider any interval (51,^2) such that B t ^ X t for all t G (si,^)- Suppose 
without loss of generality that B t < X t for all t G (si,s 2 ) Choose an arbitrarily small 
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5i > 0. Note that as 8 — > 0, the open sets W c (8) converge to the complement of {(s,x) : 
s > to,B s = x}. Let 82 > be so small that the (closed) portion of the graph of X t 
between si + 81 and S2 — 81 does not intersect Wfa)- Let sq = s\ + 8\. Since the Y t n 
converge to X t , there exists a sequence rrij such that x]^ 3 — > X So . For sufficiently large 
j, the point (so, xl^™ 3 ) lies outside VF(52) and we also have 1/rrij < 82. Then, for t in a 

^ 1 / m j 

neighborhood of so, the function x\ 3 must be given by X t 1/mj = X t s °' Xs ° . We will 
show that X t = xl°' Xa ° for t G (sq, S2 — 81). 

Suppose that this is not true and let S3 = inf{t G [sq, S2 — 8\] : X t 7^ X^°' X3 °}. 
Since (S3, X S3 ) lies outside W(52), an argument similar to the one given above shows that 
for some £3, £4 > 0, and all m > I/82, the functions X t must satisfy X t /m = X t ' 33 
for t G [S3, S3 + ^3], if \xl[ m — X S3 \ < 84. A straightforward argument now implies that 
for large n, y t n = X t 33 for t G [S3, S3 + 8s], and this in turn proves that X t = X^ 3 ' 83 
for t G [s 3 , s 3 + <y. This contradicts the definition of s 3 and proves our claim. 

Thus X t satisfies (2.3) on (si + 81, S2 — 8\) and, in view of arbitrary nature of 8\, 
the same claim extends to the whole interval (si, s 2 ). The argument applies to all intervals 
(si,S2) such that B t 7^ f{Xt) for all t G (si,S2)- This implies that X t is a Lipschitz 
solution to (2.3). The proof of the existence of a Lipschitz solution is complete. 

The existence of the solution to (2.3) for t < to may be proved in a completely 
analogous way. The two solutions can be combined into one function X t in an obvious 
way. It remains to check if the differential equation (2.3) is satisfied at t = t . It is easy 
to see that if B to < xq then dX t /dt = Fi(X t ) for all t in some intervals (to — 8, to) and 
(to, to + 8) with 8 > 0. This and the continuity of X t at t = to evidently imply that 
dXtjdt = Fi(Xt) for t = to and so (2.3) is satisfied for t = to- The case when B to > xo is 
analogous. When B to = xo then (2.3) is trivially satisfied by X t for t = to- 

Since the functions {X^ m , t > to} are adapted to the Brownian filtration = 
a(B s ,s G [to>£])) so is their lim sup, X t . It follows that the process {(B t , X t ),t > to} is 
strong Markov with respect to the filtration \Tf ' , t > to}- 

We will show that the function {Xt,t > to} constructed above is the largest of all 
Lipschitz solutions to (2.3), that is, if X t * is another Lipschitz solution, then Xt > X£ for 
all t > t - Consider any Lipschitz solution XI to (2.3) and suppose that XI > X t for some 
t > t - Then there must exist 8 = l/mj such that X? > X^ for some t > t - Fix such 8 
and let S be the infimum of those t such that X? > X^. If S G [tj + 8,tj + i) for some j, 
then Xg = Xg 7^ Bs a.s., and, by continuity, we must have X* 7^ B s and X 5 S 7^ B s for 
all s in some non-degenerate interval [S, S + 81). On this interval one of the conditions 



14 



in (2.3) is satisfied by both X* and X$, so X* = X s s = X s > x *s for all s G [S, S + 6 J or 
X* = X% = X s > x *s for all s G [S, S + <Ji). This contradicts the definition of S. Next 
suppose that S G [tj, tj + 5) for some j. On this interval, the derivative of Xf is equal to 
(3. It is is easy to see that a Lipschitz solution X^ to (2.3) cannot grow faster than that 
on this interval, and so S > tj + 5, a contradiction which completes the proof of our claim. 

A similar construction gives a solution {X t , t < to} to (2.3) which is maximal among 
all Lipschitz solutions on the interval (— oo, to] with constant (3. Note that Xt is measurable 
with respect to the a-field cr(B s , s G [t, to]) f° r t < to- 

The maximal solution Xt of (2.3) is consistent in the following sense. Consider a 
fixed path {B t , t G M} and the corresponding maximal solution X t . Now choose any s > 
and suppose that X s = z. Let {X*, u > s} be the largest Lipschitz solution with constant 
f3 for the equation (2.3) on the interval [s, oo) with the initial condition X* = z and the 
path {B t ,t G M} truncated to {B t ,t > s}. Then it is easy to see that X* = X u for all 
u > s. It follows that for s > 0, the portion {X t , t G [s, u]} of the solution to (2.3) may be 
defined only in terms of X s and {B t: t G [s, w]}. 

In a similar fashion we can construct a minimal solution to (2.3); this minimal 
solution is also adapted to the filtration of B t . Uniqueness would follow once we prove the 
maximal and minimal solutions are equal for all s a.s. □ 

Proof of Theorem 2.3. Let X + and X- be the maximal and minimal solutions to (1.1). 
By (2.3) the P x law of B t — X~(t) is mutually absolutely continuous with respect to the 
P x law of B t , so under P x , B t — X~(t) has a jointly continuous local time Lf such that 
sup x L x < oo, a.s. for each t. 



Let 



U(l) 



inf{t > : supL x > 1/(4/?)}. 



X 



If t < U(l) and a > 0, then 




Let a > and 



S = mf{t > : X+(t) - X~(t) > a}. 



15 



Since both X+ and X~ satisfy (1.1), if V = U(l) A S, 



X+(V)-X-(V) < 2(3 [ l {x -(u)<B u <x + (u))du 



Jo 



< W I l(0<B u -X-(u)<a)du 



< 2/3 / l(0<B u -X-(u)<X+(u)-X-(u))du 

Jo 




< 2a(3/(A(3) = a/2. 



Since X+(V) - X~(V) = a if £/(l) > S, we must have V = U(l). This is true for all 
a > 0, so = for t < U(l). 

Now assume that B t is strong Markov and let U(j + 1) = C/(j) + ° Qu(j)i 
j = 1,2,.. ., where 9 is the shift operator associated with the process B t . An induction 
argument using the strong Markov property at U(j) shows that X + (t) = X~(t) for t < 
U(j + 1) for j = 1,2, . . .. The continuity of B t and Lf easily implies U(j) — > oo, a.s., so 
X+(t) =X~(t) for all t > to- □ 

The proof of Theorem 2.5 will be split into several lemmas. 

For the remainder of the section, let 5 = (1 — 7)/4. Note that 5 G (0, 1/4) since 
7 G (0, 1). The constants ci, C2, . . ., in the proofs in this section may depend on 7 and 5. 

Lemma 2.14. Let a > 1, t < 1, A > 0, Ct = Jq 1(o<b s <As«) ds. Assume that condition 



(i) of Theorem 2.5 holds. There exist c\ and C2 independent of a and A such that for 
A > 0, 



Proof. First let us compute E(C t - C u \ T u ) for u G [0,t]. Let R = R(ct) = a 1/a . Note 
that R > 1, R = exp(ct _1 log a) < C3, and 



P(C t > A) < ciexp(-C2AaV(^ a+2 ' 5 )). 



1-R 



-1 



1 — exp(— log a/a) < log a/a < c^a 



-1/2 



where C3 and C4 do not depend on a as long as a > 1. 
By condition (i) of Theorem 2.5, 
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Let us examine 

"t a 

1= I , * ds. 



(s - up 

Suppose first that u < t/R. We observe, using the fact that R > 1, 

" t/R s a , /t\ a f t/R ds t a f* ds 



f S ° ds<( t Y f ds f 
Ju (s-up ~ \RJ J u (s-u)t a J u 



(s — u) 1 a J u (s — -u) 7 

fa rt-u j fa 

- %<^-t^\ (2.9) 

a J s> a 



On the other hand, in view of the inequality 1 — R 1 < c^a 1 / 2 , 

. ds 
ds<t a 

< t a 



t/R (« - U) 1 ' Jt/R (S ~ Up 

ds 



= t 



t/R (s - t/Rp 
ds 

s 7 

= C7 rt 1 - 7 (i -i?- 1 ) 1-7 

Recalling that a > 1 and combining with (2.9), 

j < CQt"" 1 " 7 

" a(i"7)/2 • 

Now suppose u>t/R. Then 

1 s a r* ds r l ~ u ds 

ds <r , a \ =t a 1 



L (S-W) 7 - J u (S-U)T Jo S 7 

= c 9 t a (t - w) 1 " 7 < c 9 t a (t - t/R) 1 ^ 

= C9 r +1 - 7 (i-i?- 1 ) 1 - 7 . 

As before, this is less than or equal to ciot a+1_7 /a (1 ~ 7 ' ) / 2 . 
Since 8 = (1 - 7)/4, t < 1 and a > 1, 

E(C t - C u | .F u ) < Cll At a+4<5 /« 25 < cnAt a+2<5 /a 5 . 

This says that almost surely the process E(Ct | JF U ) does not exceed C u by more than 
cnAt a+2S /a 5 for any u <t. In particular, 

E(C t - C T I .Ft) < cnAr +2S /a s 
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for every stopping time T bounded by t. We apply Theorem 1.6.11 of Bass (1995) to deduce 
that there exist cyi and C13 such that 

Eexp(c 12 C t a 5 /(At« +2S )) < c 13 . 

Our result easily follows from this estimate. □ 

Lemma 2.15. Given £ > 0, there exist ci, c 2 such that if a > 1, A, B > 0, B/A > £, and 
/3 = a + 5, then 

F(C t > Bt p for some t < 1/2) < a exp(-c 2 Ba s /A). 

Proof. Let t k = 2" 1 - fc /^, k = 0, 1, . . .. The process C t is increasing. So if C t > Bt 13 for 
some t < 1/2, then for some k > 1 we must have C tk _ 1 > B{t k y ■ Hence 

P(C t > St' 3 for some * < 1/2) < P(C tfc _ 1 > B(t k ) fi for some fc > 1) 

00 

<J2nCt k ^>B(t k f). (2.10) 
fe=l 

Using Lemma 2.14, this is bounded by 

c 3 exp ( - c,B{t k fa 5 /{Atttl 5 )) 

k=i 

£ c 3 exp ( - C4 ^ 2 -/ 3 - fc -(- 1 -( fc - 1 )// 3 )(«+ 25 )) 



fc=i 



00 



<5 



= J>exp (_ C4 ^2 fc ^+ 5 -(«+ 25 )/( a + 5 )). 
fc=i 

Since a > 1 and <5 G (0, 1/4), the quantity 2<M«+2<5)/(a+<5) i s bounded below and above by 
absolute constants, so the last displayed formula admits a bound 

Ba 5 



£ c 3 exp (-c 5 ^2 ^) (2.11) 

k=i 

= c 3 exp ( - c 5 ^) £ exp ( - c 5 |^(2^ - 1)). 



k=i 

The infinite sum in the last expression is bounded by 

EB, h z/a \ \ / Cz.Bk5 log 2\ 

exp(-c 5 -(2^ _ i)) < £exp ( - — ) 

1 



Af5 

fc=i fc=i 



< 



l-exp(-c 5 S51og2/(A/3)) 
< c 6 A/3/S. 
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Combining this with (2.10) and (2.11) we obtain 

Ba 



¥(C t > Bt fi for some t < 1/2) < c 3 exp ( - c 5 -^-^jc 6 A/3/ B 

= c 3 exp ^ - c 5 ^— + logc 6 + \og(A/B) + \og(a + 8) J 

( Ba s , , > , \ 

< c 3 exp ^ - c 5 ^— + logc 6 - log£ + log 2 + logaj. 



Ba 5 
Ba 5 

The last expression is less than 

( Ba s \ 
c 7 exp ^ - cg^-J 

for suitable C7 and eg (depending on £ and S) and all a > 1. □ 

Let Xf and X/~ be the maximal and minimal solutions to (2.3) constructed in in 
the proof of Theorem 2.2. Let Y t = X+ - Xf . We will show Y t = 0, a.s. for t < 1/2. 

Lemma 2.16. For each s, 

F(X+ = B a ) = 0, a.s. 
and similarly with X+ replaced by X~ . 

Proof. We know is a process whose paths are Lipschitz continuous. By assumption 
(ii) of Theorem 2.5, there exists a probability measure Q which is equivalent to P and such 
that the Q law of B s — is the same as the P law of B s . Then 

Q(X+ = B a ) = Q(B S - X+ = 0) = F(B S = 0). 

This is equal to zero by (2.4). Since P and Q are equivalent, the lemma is proved. □ 

Lemma 2.17. Y t = 0, a.s. ift < 1/2. 

Proof. The process Xf satisfies the equation 

X+ = x+ f t [F 1 (X+)l {xt>Bs) + F 2 (X+)l (xt<Bs) ] ds. 

J 

Xf satisfies a similar equation. Then, noting Lemma 2.16, 

Yt = [\Fi{X+) - F 1 (X-)]l {Bs<x -^ xt) ds 

J 

+ f[F2{Xt) -F 2 (X-)]l {x -< xt<Bs) ds 

J 

+ j\ F ^ X s) - F 2{X7)]\x 7 <B a <Xt) ds - 
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Therefore 



Yt<M^ (X+-X;)ds + 2M J l {x - <B3<xt) ds (2.12) 
= M J Y s ds + 2M J l (0<Bs _ x - <xt _ xn ds 



t rt 



s ' 

JO 



Ml Y s ds + 2M I l {0<Bs _ X7<Ys) ds. 



Recall that we have assumed that Fj is bounded by M. Hence, the process Y t is 
Lipschitz with constant 2M. Since X~ has Lipschitz paths, there exists, by assumption (ii) 
of Theorem 2.5, a probability measure Q equivalent to P such that under Q, {B s — X~ , < 
s < 1/2} has the same law as {B s , < s < 1/2} does under P. So it suffices to show that 
for any Lipschitz process Y s with constant M satisfying 

Y t <M [ Y s ds + 2M [ l {0< B s <Y s )ds, (2.13) 

Jo 

we have 

F(Y t ^ for some t < 1/2) = 0. 

Let 

D(A, a) = {Y s > As a for some s < 1/2}. 

As Y is Lipschitz with \Y t \ < 2Mt, then D(SM, 1) =0. Let e > and let rj = 1/4. We will 
choose N > 1 and jo > in a moment. Let Aj = iV J if j < j and Aj = (l + r)) J ' for j > 
jo- Let aj = l+jd. We want an estimate on the probability of D(Aj + i, ctj+i) — D(Aj, aj). 
1£w£ D(Aj, aj), then Y s < A jS a i for all s < 1/2, and so from (2.13), for t < 1/2, 

Y t < M f A jS a ids + 2M I l ( o <Ba <A jS ^) ds ( 2 - 14 ) 
Jo Jo 

MAd a > +1 n „ f' 

+ 2M / l(o<s s <^ jS ^)^- 



CXj + 1 







Let £ = (1 — 7])/2M and let ci and c 2 be constants chosen as in Lemma 2.15 
(depending on £). Find large jo so that 

(l+j 5) s / 2 /2M> 1, (2.15) 
M 
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< 77(1 + 77), (2.16) 



and 

oo 

d ^ exp(-c 2 (l - rj 2 )(l + j5) 5 ' 2 ) < e/2. (2.17) 

3=3o 

Next choose iV large so that 

AT > M/r} (2.18) 

and 

2j ci exp(-c 2 (l - rj)N/2M) < e/2. (2.19) 

For j > j , we have 

• < r}A j+1 , (2.20) 



(ay + 1) 

using (2.16). The same inequality holds for j < jo in view of (2.18). 

In view of (2.14) and (2.20), for u> to be in D(Aj + i, Oij+i) — D(Aj,aj), we must 

have, 

J* ho<B s <A jS <*i) ds > Y t /(2M) - (2-21) 

A .f(Xj + l 

>A j+1 t a '+ 1 /(2M) 1 



2( aj + 1) 
> (l- V )A j+1 t a ^/(2M) 

for some t < 1/2. Recall that we set £ = (1 — rj)/2M and note that for all j we have 
(I — 1]) Aj + i / (2M Aj) > £. By Lemma 2.15, the probability that the inequality (2.21) holds 
is less than or equal to 



Using (2.15) and (2.17) for j > jo, we obtain 

Cl £ exp ( - c.^f^a 5 ) < Cl £ exp(-c 2 ^f (1 + rf /2 (l + ^) 5/2 ) < e/2. 



From (2.19), 

^ (1-^)^7+1 A ^ / (\-ri)N\ 

3=0 3 3=0 

< 2j oCl exp(-c 2 (l - r])N/2M) < e/2. 

Hence, 



3=0 ~ " 3 
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and so 

oo 

F((J D(A jtaj j) <e. 

3=0 

If u i \Jf = oD(Aj,OLj), then Y t (u) < A 3 t^ < (1 + v) j N jo (1/2) 1+ ^ for all j > j 
and all £ < 1/2. Since (1 + 77) (1/2) < 1, letting j — > 00 shows Yt(u;) = 0. Therefore 

P(lt ^ for some £ < 1/2) < e. 

Since £ is arbitrary, this proves the lemma. □ 



Proof of Theorem 2.5. By Lemma 2.17 we have Y t = a.s. for £ < 1/2. If we consider 
the law of B t+1 / 2 given .F1/2, it is not hard to see that assumptions (i) and (ii) of Theorem 
2.5 apply to this process as well. So we apply the same argument to X^ +1 ^ 2 and X~ +1 ^ 2 , 
and we obtain Y t+1 / 2 = for t < 1/2, or Y t = for t < 2(1/2). By an induction argument, 
we then have Y t = for all £, which proves uniqueness. □ 



Proof of Theorem 2.6. The existence and strong uniqueness of solutions X* to (2.6) 
can be proved in the same way as in Theorem 2.1. 

Consider any sequence e n [ and with a slight abuse of notation let X™ = X\ n . 
Since all functions t — > X™ are Lipschitz with constant /?, we may suppose, passing to a 
subsequence, if necessary, that X™ converge to a function X^° . In order to finish the proof, 
it will suffice to show that X^° = X t . Since the equation (2.5) has a unique solution a.s., 
it will be enough to show that if uj is not in the null set where uniqueness does not hold, 
then X^°(oj) is a solution to (2.5). The functions X™ are Lipschitz with constant /?, so the 
same is true of X^°. Let A be the set of times t such that X^° = B t . The complement of 
the set A consists of a countable number of open intervals. Let I = (ti,t 2 ) be one of the 
intervals in the complement of A. Fix any £3 e / and suppose without loss of generality 
that X^ > B t3 . Choose some £4 G (£1,^3) and £5 G (£3,^2) and let a be the infimum of 
X t °° — B t over t G (t 4 , £5). For sufficiently large n, we have £ n < a/3 and |X t n — X t °°| < a/3 
for all t G (£4, £5). It follows that for large n and t G (£4, £5), we have X™ — B t > a/3 > e n . 
Hence, for such n and £, a'X t n /a't = (3 2 . This shows that dX^/dt = (3 2 for all £ G /. The 
same argument works for all other intervals in the complement of A. There is nothing to 
check for t G A, so X^° is a solution to (2.5). □ 
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3. Local time. In the remaining part of the article we assume that B t is a Brownian 
motion. In this section we will exclusively deal with solutions to (1.1). We will find several 
explicit formulae for the local time spent by B t on the paths of the process X t . Moreover, 
we will prove analogues of the Trotter and Ray-Knight theorems. The results on local 
times provide information about the behavior of the function y — > Xf , for fixed t; see 
Remark 3.9. 

The first part of the section deals with exit systems. Some of our results on exit 
systems may be of independent interest. We refer the reader to Blumenthal (1992), Burdzy 
(1987), Maisonneuve (1975) or Sharpe (1989) concerning the fundamentals of excursion 
theory. 

In this section, we will assume that to = and study the portion of the solution X t 
to (1.1) for t > only. 

Let D = {(&, x) G M 2 : b = x}. We will construct an exit system (H x ,dL) for the 
process of excursions of (B t ,X t ) from the set D. The first element of an exit system is a 
family of excursion laws H x . An excursion law H x is an infinite cr-finite measure on the 
space C* of functions (ej, e 2 ) defined on (0, oo) (note that is excluded) which take values 
in R 2 U {A}. Here A is the coffin (absorbing) state. Let v be the lifetime of an excursion, 
i.e., v = inf{t > : (e t \e 2 ) = A}. Then fP-a.e., we have (ej,e 2 ) G M 2 for t G (0, u) and 
(ej, e 2 ) = A for t G [v, oo). The measure H x is strong Markov with respect to the transition 
probabilities of the process {(B t , X t ),t > 0} killed at the hitting time of D. Moreover, 
the IP-measure of the set of paths for which limtj.o(eJ, e 2 ) ^ (x,x) is equal to 0. The 
second element of the exit system, dL, denotes the measure defined by a non-decreasing 
process L t . The process L t is a continuous additive functional, also known as a local 
time, for (Bt,Xt) on D. The process Lt does not increase on any interval (s,u) such that 
(Bt, Xt) ^ D for t G (s, u); that is, L s = L u for such intervals. Consider a maximal interval 
(s,u) such that B t 7^ X t for t G (s,u). Suppose L s = r. Let (ej,e 2 ) r = (B s+t , X s+t ) for 
t G (0, u — s) and (ej, e 2 ) r = A for t > u — s. Let fi(r) = mf{t > : L t = r}. The collection 
of all "excursions" {(r, (e 1 ,e 2 ) r )} is a Poisson point process which, roughly speaking, has 
random mean measure (r2 — r±) H^ r \A)dr on the set (ri, r^) x ^4- 

Next we apply some transformations to the excursions and excursion laws in order to 
simplify our description of the exit system. First, we note that by the translation invariance 
of the Brownian motion B t and equation (1.1), the distribution of (e\ — x,e 2 — x) under 
H x is the same for every x G KL Let this distribution be called Hi. For i^i-almost all 
excursions, the second component e 2 is a linear function of t until the excursion lifetime 
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v, with the slope equal to f3\ or (3 2 . In the first case, e\ > e\ for t G (0, v), while the 
inequality goes the other way in the second case. Let H 1+ denote the part of the measure 
Hi which is supported on excursions with e\ > ef and let Hi- be the part supported on 
the set where el < ef . Let H 2+ be the distribution of {el — ef,t G (0, z/)} under Hi + and 
let H 2 - have the same definition relative to Hi-. Note that, by definition, the excursion 
laws H 2+ and H 2 - are supported on paths in MU {A} rather than M 2 U {A}, since the 
second component becomes irrelevant after our last transformation. 

Our transformations preserve the strong Markov property, but the last transforma- 
tion creates a drift so that the measure H 2+ has the transition probabilities of Brownian 
motion with drift — fli, killed upon hitting 0. It is standard to show (see, e.g., Theorem 
4.1 of Burdzy (1987)) that for any event A defined in terms of the process after some fixed 
time so > 0, we have, up to a multiplicative constant, 

H 2+ {A) = \ir*J-Q*_ (}i {A), (3.1) 

where Q-^ stands for the distribution of Brownian motion with drift — Pi, killed at the 
hitting time of 0. The normalization of the excursion laws is arbitrary as long as it matches 
the normalization of the local time, so we can use the normalization in (3.1). We next 
choose the normalization of the local time so that it matches that of H 2+ . Given the 
normalization for H 2+ , the normalization for H 2 - is no longer arbitrary and we will have 
to prove that 

H 2 -(A)=lim^-Q^ 2 (A). (3.2) 

Unless specified otherwise, all excursion laws in this paper will be normalized as in 
(3.1) or (3.2). 

Let H 3 denote the excursion law for excursions of Brownian motion without drift 
away from 0. Let us split H 3 into positive and negative parts H 3+ and H 3 - : as in the case 
of H 2 . We normalize H 3 using a formula analogous to (3.1). Recall that v denotes the 
lifetime of an excursion e, and that (3.1) defines the normalization of H 2+ . 

Lemma 3.1. (i) On the set where v < oo, 

dH 2+ 



dH 



3+ 



(e) = exp(-/3 1 ^/(2 ( x 2 )). 



(ii) For a fixed time s G (0, oo), the conditional distributions of H 2+ and H 3+ given 
{u = s} are identical. 
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(Hi) If fa < then H 2+ (u = oo) = 21/^l/a 2 . 

(iv) Formula (3.2) is the correct normalization for H 2 - 



Parts (i)-(iii) of Lemma 3.1 have obvious analogues for H 2 _. 

Proof. Fix arbitrary < so < s± < oo and let A be an event measurable with respect to 
cr{et,t G (sojSi)}. Since Hs + is assumed to be normalized using a formula analogous to 
(3.1), we have 

H 2+ (A n{v = si}) _ Hm Q- Pl (A n {u = si}) 

H 3+ (An{v = Sl }) ™ Q*(An{v = Sl }) 

An application of Girsanov's Theorem, as in Karatzas and Shreve (1988) ((5.11), p. 196), 
shows that 

O x n (An {v = sA) 

= exp(^ 1 /a 2 -/3 1 2 s 1 /(2a 2 )). 



exp(-/3 1 2 Sl /(2 ( r 2 )), 



QUAn{u = Sl }) 
This and the previous formula imply 

H 2+ (An{u = Sl }) 

H 3+ (An{is = Sl }) 

which then easily implies (i) and (ii). 

As for (iii), we start with the formula 

Qlfriv = oo) = 1 - exp(2x/3 1 /a 2 ), 

with /?i < (Karlin and Taylor (1975), p. 362). Then (3.1) yields 

H 2+ (u = oo) = \im^Q x _ Pi (v = oo) = 2\p t \/a 2 , 

xlO \X\ 

as desired. 

It remains to prove (iv). Fix arbitrarily small 7 > and let 
Ax = A x (t) = {max \X S \ > t 1/2+ ~<}. 

s<t 

Note that \X t \ < f3t < ^2+7 for small t > so we have P(Ai(t)) = if t is small. However, 
we will prove the result using only the property that \im t ^ ¥(A 1 (t)) = because we will 
need this version of the proof later in the paper. Let us take to = and xq = so that 
Xq = 0. Note that the excursion law normalization does not depend on to and xq. Let 
A + = A + (s) be the event that the first excursion (ej,e 2 ) of (B t ,X t ) from D with the 
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property that \e\ — e 2 | > s 1 / 2 " 1 " 7 / 2 for some t G (0, z/), also has the property that e\ > e 2 
for t E (0, v). Let A_ be the analogous event with e\ < e 2 . Let T(a) be the hitting time 
of a by For small s > 0, 

{T(s l/2 +7 /2 + s l/2 +7) < T( _ s l/2 +7 /2 + s l/2+ 7) < s} c A+(s) y Ai(s); 

and 

{T (_ s l/2+7/2 _ s l/2+ 7) < T(s l/2 +7 /2 _ s l/2+ 7) < g} Q A _ {g) y 

It is elementary to check that 

lim P(T( S 1 / 2 + 7 / 2 + 5^2+7) < T (_ s l/2+ 7 /2 + s l/2+ 7) < s) 
s— >0 

= lim P(T(-S 1 /2+7/2 _ s l/2+ 7) < T(s l/2 +7 /2 _ s l/2+ 7) < g) = 1/2 _ 
s— >0 

This and the fact that lim t ^ P(Ai(t)) = imply that 

limP(A+(s)) = limP(A_(s)) = 1/2. (3.3) 

s^O s^O 

The scale function S(y) for Brownian motion with drift —f3\ is given by 
S(y) = exp(2/3iy/a 2 ) (Karlin and Taylor (1981) Chapter 15.4). Let F h be the event 
that the difference between the maximum and the minimum of an excursion exceeds h. 
Then, by (3.1), 

H 2+ (F h ) = lim -Q^ (T h < T ) = lim - • f X J - M 
xio x M1 x[Q x S[h) — S(Q) 

1 exp(2/Wa 2 ) - 1 2/?i 1 

= lim 



xio x ex.p(2/3 1 h/a 2 ) - 1 a 2 ex.p(2/3 1 h/a 2 ) - 1 ' 

An analogous formula holds for H2-(Fh), but we will write it with an additional multi- 
plicative constant c\, since we have not proved that (3.2) is the right normalization yet: 

H 2 .(F h ) = Cl 2f32 1 



a 2 exp(2/3 2 /i/^ 2 ) - 1' 

Our goal is to show that c\ = 1 is the correct choice for the constant. 

Excursion theory tells us that the arrival times for excursions (e^,e 2 ) of (B t ,X t ) 
from D with the property that \e\ — e 2 | > s l ^ 2+1 ^ 2 for some t e (0, u), and with e\ > 
ef for t E (0, v), form a Poisson point process on the local time scale with intensity 
H2+(F s i/2+t/2). This process is independent from the analogous process of excursions 
with e\ < e 2 . Formula (3.3) tells us that for small s, the probability that the first arrival 
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for the first process is earlier than the first arrival for the second process is close to 1/2. 
Hence, the ratio of the intensities for the two Poisson point processes must converge to 1 
as s — > 0. Therefore, we must have 



Remark 3.2. (i) Lemma 3.1 (iv) can be used to prove uniqueness for (1.1). In order to 
do so, one would have to consider an exit system for the maximal Lipschitz solution X t 
to (1.1), constructed as in the proof of Theorem 2.2, and the analogous exit system for 
the minimal Lipschitz solution. Lemma 3.1 (iv) shows that both exit systems are identical 
but this can be true only if the maximal and minimal solutions are the same. We will not 
formalize this argument as it cannot be easily generalized to non-Markov processes. The 
delicate part of the argument would be to show that the maximal solution X t is the sum 
of the excursions and it contains no component corresponding to a "push" proportional to 
local time. 

(ii) According to Lemma 3.1 (i)-(iii), if f3\ < 0, the excursion laws H%*. and 
agree on the set of excursions with finite lifetime and the only difference is that H^X gives 
some mass to excursions with infinite lifetime, while 1 does not. 

For every x G R consider the solution Xf to (1.1) with Xq = x. Let denote 
the local time of Yf = B t — Xf at 0, defined earlier in this section as the local time of 
(B t ,X£) on the diagonal, accumulated between times and t. Note that this is not the 
local time of a one-dimensional diffusion at level x. 

Proposition 3.3 If (3\ > and f3 2 < then for every 



Proof. Fix some x G R. Our assumptions that (3\ > and f3 2 < imply that there will 
never be an excursion of B t from Xf with infinite lifetime, since the drift will always push 
the excursions of Yf towards 0. This in turn implies that will grow to infinity a.s. 




However, this is possible only if ci = 1. This completes the proof of (iv). 



□ 




a.s. 
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Recall that we used Q v _p x to denote the distribution of Brownian motion with drift 
—Pi, killed at the hitting time of 0. By Theorem 7.5.3 of Karlin and Taylor (1975), we 
have for y > 0, 

Q v s(ue dt) = — !^L= exp ( J\v\-W) 2 ) ^ 



where v denotes the lifetime of the process. The same formula holds for y < 0, with Pi 
replaced by Pi- Using (3.1)-(3.2) we obtain 

1 1 

B 2 {y E dt) = lim — Q y R (v e dt) + lim --Q y * (u e dt) 

= lim - — == exp _ at + hm - — == exp _ dt 

= exp[-(tf /2^)t]* + eM-(Pl/^)t]dt. 

Let V/ be the inverse local time, i.e., Vf = inf{t > : Lf > s}. The process Vf is the 
sum of lifetimes of excursions which start before the local time reaches the level s. The 
Poisson character of the excursion process easily implies that 

poo 

EV S X = s tH 2 (vE dt) 
Jo 

= S / t ,3/2 fn~ exp[-(ft/*r 2 )t]dt + s / t exp[-(p 2 2 /2a 2 )t]dt 
Jo (rt 6 i' i y/2TT Jo at^/^yz-K 



' +4t)«- ( 3 - 4 ) 



This and the strong law of large numbers for the Levy process s — > Vf (see p. 92 of Bertoin 
(1996)) imply that 



a.s., as s — > oo. This can be easily translated to the statement of the proposition. □ 



We note that if Pi, P 2 > 0, then we will eventually have > B t , for every x. 
Hence, in this case, < oo for every a.s. We will prove the next lemma under the 

assumption that Pi — P 2 > 0. We believe that similar statements hold when Pi — P 2 < 
but technical difficulties prevent us from giving a formal proof in that case. 

The following lemma contains the most complicated and technical argument in the 
whole article. 
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Lemma 3.4. (i) Fix x, a, Pi, p 2 > and assume that (3i — p 2 > 0. Then 

E(L^ S \L% = o) = a- S-^-(l - exp(-2a(/3 1 - (3 2 )/a 2 )) + o(S), 

Pi — P2 

for 5 | 0. 

(ii) If a, p u p 2 > 0, x < 0, and pi - p 2 > then 



= a) = a + 5 



fa 



fa 



fa- fa fa- fa 



exp(-2a(/3 1 -/? 2 )/(7 2 ) 



for 8 t 0. 



Proof, (i) Recall that So = 0, that we have fixed x,a,fa,fa > and assumed that 
fa — fa > 0. Since the proof of the lemma is quite long, we will split it into several steps. 

Step 1. We start with some transformations of the processes Xf and B t which will enable 
us to look at from a slightly different perspective. It is perhaps not necessary to make 
these transformations, but we find the transformed problem much easier to comprehend 
than the original one from an intuitive point of view. 

We first offer a rough guide to our notation (whose validity is limited to this proof). 
Different Brownian motions with different drifts and reflected barriers will be denoted B\ , 
for j = 1,2,.. .. The notation L\~ and L J t will refer to the local time of Bf on the lower 
and upper reflected barriers (if any). We will write Vj(a) = inf{t : L\~ = a}. 

It is well known that the set {t : B t = 0} has zero Lebesgue measure, so the 
Girsanov theorem implies that the same is true of the set {t : Xf = Bt}. We will excise all 
intervals where Xf > B t . First, we define a clock Ci(t) = <B s }ds and its inverse 

bi(t) = inf{s : Ci(s) > t}. Since fa, fa > 0, we will eventually have Xf > B t , so we let 
ui = sup{Ci(t) : t > 0}. Then we define new processes on the random interval [0,ui] by 

B bl{t) -faib^t) -t), 

xz i{t) -fa(h(t)-t), 
x^ { 5 t) -fa(h(t)-t). 

For t G [0, u\], we have X\ ,x = x + fat, the process B\ is a Brownian motion staying above 
and reflected on the line t — > x + fat, and the process x] ,x+S is a solution to (1.1) with 
B t replaced by B\. 



X 



X 



l,x+(5 
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Next we similarly excise the intervals where xl' x+S < Bf. Let us define a new 
clock C 2 (t) = J 1 inverse 62 (t) = inf{s : C2(s) > t}, a random time 

U2 = sup{C 2 (t) : t > 0}, and processes 

5? = Bj, (t) - ft (fc (*)-*). 
^ = ^-^(62(0-0, 

^ 2 ' X+5 =^( t | 5 -/3i(^(t)-t). 

For t G [0, U2], we have X^' x = x + Pit and xf ,x+S = x + 5 + fat. The process {B%, t G 
[0, W2]} is a Brownian motion reflected on the lines t — > x + /3it and t — > x + 5 + and 
confined to the region between them. Note that Bq = x a.s. and that the lines t — > x + 
and t — > x + 5 + /^t intersect at t = 5/ (Pi — P2) so necessarily u 2 < 5/ — /3 2 ). 

The time u 2 corresponds to the start of the infinite excursion of B t below the graph 
of Xf . By excursion theory and Lemma 3.1 (iii), the distribution of L^ 2 is exponential with 
mean cr 2 /(2/3 2 ). Hence, we may assume that the process Bf is generated in the following 
way. Suppose that Bf is a Brownian motion starting from Bq = x, reflected on the lines 
t — > x + I3\t and t — > x + 5 + fat and confined to the region between them, but defined for 
all t G [0, 5/ (Pi — P2)) rather than confined to some random time interval. Let Lf~ be the 
local time of Bf on the line t — > x + Pit and let Z be an exponential random variable with 
mean cr 2 /(2/3 2 ), independent of Bf. If Vs(s) = inf{t : Lf~ = s}, then the distributions of 
the processes {Bf,t G [0,-u 2 ]} and {Bf,t G [0,i> 3 (Z)]} are the same. 

Let Lf + be the local time of Bf accumulated on the line t — > x + 5 + /3 2 t. The 
distribution of L^ 5 given {L^ = a} is the same as the distribution of L^ a y so we will 
try to find an approximate formula for EL^^ . 

We continue our transformations. Let Bf = Bf — x — Pit. The process Bf is a 
Brownian motion starting from 0, with drift — Pi, reflected on the horizontal axis and the 
line t — > S — (Pi — P 2 )t. The processes L\~ and Lf + can be identified with the local times 
L\~ and Lf + of Bf on the horizontal axis and the line t — > S — (Pi — P2Y, resp. Hence, it 
will suffice to show that the estimate given in part (i) of the lemma holds for EL^ a y 

Step 2. In this step we will obtain some estimates for reflected Brownian motions using 
excursion theory. Let Bf be a Brownian motion with drift — Pi , confined to positive values 
by reflection on the horizontal axis. The Green function G(z, y) for Brownian motion with 
drift —Pi, killed upon hitting is given by 



G(z,y) = ±- 



1 2 ^ z 1 1 
exp ( J - 1 
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exp 



2Piy 

a 2 



for < z < y < oo, by (3.15) in Section 15.3 and Section 15. 4. B of Karlin and Taylor 
(1981). Let G 5 H (y) denote the Green function for the excursion law H 5 of B\ from 0, i.e., 
the function defined by 



l{e(t)e[z 1 ,z 2 ]}dt 



y-s: 



A formula analogous to (3.1) yields 



G 5 H (y) = \un-G(z,y) 

zj.0 Z 



exp 



G 5 H (y)dy. 



for y > 0. 

Consider some 5i > and excise excursions of above the level S\ , just as we did 
with the excursions of B t and B\. Let C 3 (t) = J* l{ B s< Sl yds, 6 3 (t) = inf{s : C 3 (s) > t}, 
and B$ = B^ t y The process Bf is a reflected Brownian motion in [0, <5i]. Let G%(y) 
be the Green function for the excursion law H 6 of Bf from 0. It is clear from the nature 
of the transformation which generates Bf from the paths of £?f that G%(y) = G 5 H (y) for 
y G (0, 8i). Hence, 

r*i r*i o / 2/3^ \ . 1 



ff fl (i/)= I' G%(y)dy= /'^exp 



1 — exp 



Let and L® + denote the local time of Bf at and Si, resp. Let vq(s) = mf{t : L\~ = s}. 
The random variable vq(s) is the sum of the lifetimes of excursions of Bf from which 
occur before L\~ reaches the level s. The last formula and excursion theory give 



Kvq(s) = s 



1 

Pi 



exp 



2(h6i 



= ST](5i). 



(3.5) 



Next we will derive an estimate for H§{v > t). Recall that Q z _p 1 denotes the 
distribution of Brownian motion with drift —fii, killed upon hitting 0. Let Q-p 1 denote 
the distribution of Brownian motion starting from z G (0, <5i), with drift — j3\, reflected at 
8i, and killed upon hitting 0. It is easy to see that 

Q z _p 1 {v>t)<Q z _ /3l {v>t), 

for alH > and z G (O,^). By Lemma 3.1 (i), Theorem 5.1 (iii) of Burdzy (1987), and 
scaling, 

H 6 (u >t)< lim -Q Z _ B , iy > t) 

zlO Z 11 



1 f°° 1 

< lim -Q Z _ B (v>t) = H 5 (v>t)< / -(2irs 3 )- 1/2 ds. 

' zio z Pl J t a 



(3.6) 
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A simple argument based on scaling and the Markov property applied at times t = k5j, 
k = 1, 2, . . ., shows that there exists a constant c\ > 0, such that 

Q-pM > *) < exp(- Cl ta 2 /^), (3.7) 

for all t > Sf/a 2 and z G (0,S\). Another standard estimate is 

Q- Pl ^>Sl/a 2 )<zc 2 . 

This combined with the previous estimate gives (with possibly new values for the con- 
stants) , 

Q z -p^{y > t) < zc 2 ex.p(-cita 2 / 8l) , 

for all t > Sf/a 2 and z G (0,6\). We obtain from this an estimate analogous to (3.6) but 
applicable for t > Sf/a 2 : 

H 6 (u >t) <c 2 exp(-cit(T 2 /^). (3.8) 
Since the excursion process is a Poisson point process, we have from (3.6) and (3.8), 

POO 

Varv 6 (s) = s t 2 H 6 (v G dt) 
Jo 

< s f 1 ^ -t 2 (2nt s )- 1 / 2 dt + s -(5 2 /a 2 ) 2 (2nt s )- 1 / 2 dt 

J0 a Jsf/a 2 & 

f°° 5 2 

+ s t 2 c 2 ^-^ex.p(-c 1 ta 2 /5j)dt 

J 5 2 /(T 2 C X O l 

< c 3 s5l/a 4 + c 4 s5f/a 4 + c 5 s5f/a 8 < c 6 sSf. (3.9) 

Step 3. We will find a link between processes reflected on sloped lines (in space-time) and 
within an interval. We will need to define some more variables. First of all, sq > should 
be considered a small constant whose value will be chosen later in the proof and which 
does not change with 5. Recall r\ defined in (3.5). Let uq > and Si G (0, 5) be defined by 
the following two equations uq = (S — 5i)/ (fli — (3 2 ), and sq = uo/rj(Pi, 5, a 2 ). 

Recall that Bf is a reflected Brownian motion in [0, 5{\ and note that now 5i is 
defined relative to 5. Let Bj be the analogous reflected Brownian motion in [0,6]. 

Note that 5 — (fix — f3 2 )t > 5i for t G (0, u ). Hence, on the interval (0,-u ), the 
upper reflecting boundary for Bf lies below that for Bf. This relationship between the 
upper reflecting boundaries implies that the excursion measure distribution of the lifetime 
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of an excursion from of the process Bf is stochastically larger than that for an excursion 
of Bf, for excursions within the interval (0,-uo)- It follows that one can construct Bf and 
Bf on a common probability space so that ve(s) A t < i>4(s) A t for all t < uq, where 
v±(s) = inf {t : L\~ = s}. On the other hand, 5 — (Pi— f3 2 )t < S for t > 0, so the analogous 
relationship for B 7 goes in the opposite way, i.e., v 7 (s) A t > i><t(s) A t for all t > 0. 

Although the process Bf starts from 0, by construction, it will be necessary to 
consider the case when it starts from some other value; in other words, we will now 
consider a process with the same transition probabilities but a different starting point. 
The starting point y will be reflected in the notation by writing ¥ y or E y , as usual. 

Let Tq be the hitting time of for the process B\ for % = 4, 6, 7. By the previous 
remarks, we can construct versions of Bf and Bj on the same probability space so that 
they start from the same point y and Tq A t < Tj A t for t < uq- 

By the strong Markov property applied at Tj , we have Ey V7 (s) = E^T 7 + E°v 7 (s). 
It follows easily from (3.7), applied to 5 rather than Si, that 

E v Tq < c 7 5 2 /a 2 . (3.10) 

Consider arbitrarily small e G (0, 1/4). We obtain using (3.5) (applied with 5i replaced by 
5) and (3.10), 

E> 4 (so)) < ^ y (v 7 (s )) < E y T 7 + E°(v 7 (s )) < c 7 5 2 /a 2 + u . 

For small 5 > 0, (3.5) shows that t](d) is approximately 26 jo 2 . Hence, uq = sorj(S) is 
approximately equal to 2sq5/ct 2 . This shows that for small 5, the last displayed inequality 
yields 

E y (v 4 (s )) < uo(l + e) = s v(S)(l +e)< 2s S(l + e) 2 /a 2 . (3.11) 
Next we will find a lower bound for the same quantity. 

By the strong Markov property applied at T 6 , we have E y v 6 (s) = E v Tq + E°v 6 (s) 
and Var (v 6 (s) \ B$ = y) = Var (T 6 | B% = y) + Var (v 6 (s) \ B$ = 0). We have an estimate 
analogous to (3.10): 

E v Tq < c 7 5i/a 2 , (3.12) 
and another estimate following from (3.7): 

Var(T 6 |B 6 = ?/ )< C8 ^/a 4 , (3.13) 

for any y G [0,5i]. 
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We obtain using (3.5) and (3.12), 
W(v 6 (s (l - £))) < JE^Tq + E°(v 6 (s (l - e))) < c 7 5j/a 2 + s (l - e)^). 

For small 5 > 0, 5i is also small and (3.5) shows that 77(^1) is about 25\/a 2 . Hence, 
so(l — e)t](5i) is approximately equal to 2so(l — e)8i/a 2 . This shows that for small 5, the 
last displayed inequality yields 

W(v 6 (s (l - e))) < s (l - 6/2)71(6!) < s (l - e/2)7/(<y) = u (l - e/2). (3.14) 

A similar estimate for the variance follows from (3.9) and (3.13), for small 5, 

Var (v 6 (s (l - e)) \ B$ = y) < c 8 5 4 /a 4 + c 6 s (l - e)5\/a 4 < c 9 s (l - e)8f/a 4 . (3.15) 

This estimate, (3.14) and the Chebyshev inequality yield, 

P"M*o(l - e)) > no) < c ^~e)5f/a 4 ^ c 10 5? c 10 5? 



(euo/2) 2 ~ ri(5)e 2 u o 4 r] 2 (5)e 2 s a 4 ' 
For small 5 we have 

5/a 2 < 7](6) < 45/a 2 . (3.16) 

Hence, 

PWl- £ ))> U o)< w J;4 og4 <^. (3.17) 
We have from (3.14)-(3.17), for small 5, 

E> 6 (s (l - £))) 2 = Var (v 6 (s (l - e)) \ B 6 = y) + (Wv 6 (s (l ~ e))) 2 
< c 9 s (l - e)5f/a 4 + (s (l - e/2) V (S)) 2 
<32s 2 (l-e/2) 2 S 2 /a 4 . 
We use this estimate, (3.5) and (3.16)-(3.17) to obtain, for sufficiently small 5, 

E>4(s )) >W(v A (s ) A«o) 
>E y (v 6 (s ) Au ) 
>Ey(v 6 (s (l-e))Au ) 

> Ey(v 6 (s (l - e))) - E^ [v 6 (s (l - e))l { „ 6 ( flo( i-e))>«o}] 

> 8 (l-e) V (Si) - (E»(«8( ao (l -£))) 2 E y (l W(so(1 _ £)) > Mo} ) 2 ) 1/2 
= So (l - e)»7(*i) - (E»(i*(a (l - e))) 2 Vy(v G (s Q (l - e)) > u Q )) 1/2 

> 2s (l - e) 2 b x /a 2 - ( 32s 2 (l - e/2) 2 5 2 /a 4 ■ ^) . (3.18) 
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It follows from (3.16) and the definition of 5i and u that for small 5 > 0, 



5 1 = 5 - s v{6)Wi - AO > 5(1 - 4s (/?i - AO/* 2 )- 

We will choose sufficiently small s > (relative to a, Pi,P 2 and e) so that Si > 6(1 —e/2). 
Then the last inequality and (3.18) yield for small 8, 

E y (v 4 (s )) > 2s (l - e) 3 V^ 2 - (3-19) 



Step 4- We will apply induction in order to obtain estimates for EPv^jso) with integer 
j > 1. At the time i>4(so), the distance between the reflecting barriers for Bf is equal to 
6 = 6 — (Pi — f32)v4(so), which is less than 5, so we can use the estimates (3.11) and (3.19) 
with 6 replaced by 5, assuming that 5 itself is sufficiently small for the estimates to hold. 
By the strong Markov property, 

E°(v 4 (2s ) - v 4 (s ) | v 4 (s )) < 2Ss Q (l + ef/a 2 = 2[6 - (pi - p2)v 4 (s )]s (l + ef/a 2 , 



and so 

¥?(v 4 (2s ) - v 4 (s )) < E°2[5 - (Pi - P 2 )v 4 (so)}s (l + ef/a 2 

< 2[S - (Pi - P 2 )2Ss (l + e) 2 /a 2 ]s (l + e) 2 /a 2 
= 2S(s /a 2 )[(l + e) 2 - 2(pi - P 2 )(s /a 2 )(l + e)% 

It follows that for small 6 > 0, 

EV(2s ) = E°v 4 (s ) + E ^ 4 (2 S o) - v 4 (s )) 

< 2S(s /a 2 )(l + e) 2 + 2S(s Q /a 2 )[(l + e) 2 - 2(pi - p 2 )(s /a 2 )(l + e) 4 ] 
= 2S(s /a 2 )[2(l + e) 2 - 2(pi - p 2 )(s /a 2 )(l + e)% 

More generally, 

E°v 4 ((j + l) So ) = E°v 4 (js ) + E°(v 4 ((j + l)s ) - v 4 (js )) 

< E°v 4 (js ) +2[5- (Pi - p 2 )E°v 4 (jso)]s (l + e) 2 /a 2 

= E Q v 4 (js )[l - 2(pi - p 2 )(s /a 2 )(l + e) 2 ] + 2S(s /a 2 )(l + e) 2 . 
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From this we obtain by induction, 

EV(jso) < E%(s )[l - 2(A - p 2 )(s /a 2 )(l + e) 2 ]^ 1 

+ 25(s /a 2 )(l + e) 2 £[1 - 2(ft - [3 2 )(s /a 2 )(l + e) 2 ] k 

k=0 

< 25(s /a 2 )(l + e) 2 [l - 2{f3 1 - (3 2 )(s /a 2 )(l + e) 2 ]^ 1 

3-2 

+ 25(s /a 2 )(l + ef ]T[1 - 2(ft - [3 2 )(s / V 2 )(l + ef] k 

k=0 

OM, /^Wi i ^2 l-[l-2(/3i-/3 2 )(Va 2 )(l + £ ) 2 p 
= 2 ' (Va Kl + £) l-[l-2( A -ft)(^)(HeF] 

= Tr^-a - [i - 2(A - /3 2 )(V^ 2 )(i + ^) 2 P)- 

Pi — P2 

Now fix an arbitrary a > 0, an arbitrarily small e > 0, and choose a sufficiently 
small small so > so that 8\ > 5(1 — e/2), and such that for some integer j we have 
jso = a, and, moreover, j is sufficiently large to imply the following: 

-^-(1 - [1 - 2{f3 1 - P 2 )(s /a 2 )(l + e) 2 Y) 

Pi — P2 

d (l-[l_2(/3 1 -/3 2 )(V^ 2 )(l + £) 2 ] a/so ) 



< " exp(-2a(/3 1 - + e) 3 /^))- 

Pi — p 2 

Then for sufficiently small 5 > we have, 

E%(a) = E%(js ) < - r ^- r (l-exp(-2a(/3 1 - [3 2 )(1 + e) 3 /o 2 )). (3.20) 

Pi — P2 

A completely analogous argument using (3.19) in place of (3.11) yields 

E%(a) > -^-(1 - exp(-2a(A - (3 2 )(1 - e) 4 /a 2 )). (3.21) 

Pi — P2 

Since e > is arbitrarily small, a standard argument based on (3.20)-(3.21) gives for 8 I 0, 
E V(a) = TT^-rt 1 - exp(-2a(A - /3 2 )/a 2 )) + o(6). (3.22) 

Pi — P2 

Step 5. The last part of the proof exploits a relationship between local time and certain 
stopping times. Recall the local times L\~ and L^ + , introduced earlier in the proof. We 
have for some standard Brownian motion Bf , 

Bt = B*-p lt + Lt-L$ + . 
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One has to check that the normalization of the local time, defined relative to the normal- 
ization of the excursion laws in (3.1), is the correct one for the above "Levy formula." This 
can be done, for example, by comparing our normalizations with those in Theorems 3.6.17 
and 6.2.23 in Karatzas and Shreve (1988). 

Note that the cr-fields generated by B 4 and Bf are identical so v±(a) is a stopping 
time for Bf . We have B 4 , x = and L 4 ~, , = a, so 

t V4(a) V4(a) ' 

= < (a) - p lV4 (a) + a - L 4 v + a) . (3.23) 

Since v±(a) is bounded by 8/((3\ — fe) the optional stopping theorem yields EB^^ = 0, 
and so, using (3.22), 

EL 4 + (a) = a - Ep lV4 (a) =a- ^7^(1 - exp(-2a(/?! - /3 2 )/a 2 )) + o(8). (3.24) 

Now recall that L 4 ^ a ^ has the same distribution as L^ 5 given {L^ = a}. This observation 
and the last formula complete the proof of part (i) of the lemma. 

(ii) The proof of part (ii) of the lemma uses a formula analogous to (3.24), but 
requires some additional work. 

Recall that 8 was positive in part (i) of the proof; it will be negative in the present 

part. 

Recall the transformations of B t from the proof of (i). It is easy to see that analogous 
transformations in the current case do not lead to B 4 which is a Brownian motion starting 
from 0, with drift — reflected on the horizontal axis and the line t — > 8 — {(3\ — $2)t 
(with 8 > 0), but instead they give a Brownian motion B 4 starting from 0, with drift — /?2, 
reflected on the horizontal axis and the line t — > 8 + (/3i — (i^t (with 8 < 0). 

A subtle but significant difference from (i) is that the infinite excursion of Bt from 
the graph of Xf will go in the direction of the graph of X® +s and so it will generate 
some more local time. By Lemma 3.1 (i) and Remark 3.2 (ii), the excursions with finite 
lifetimes have the same intensities for Brownian motions with drifts /?2 and —02, so we 
can use estimate (3.24) for the portion of the local time generated before the last, infinite 
excursion of B t from the graph of Xf . The estimate has to be modified as j3\ has to be 
replaced by /?2, and so we obtain 

« " l 5 l7^V(l - eM-MPi ~ /3 2 )/a 2 )) + o(8). (3.25) 

Pi — P2 

To this we will have to add the local time spent by B t on the graph of X^ +s during its 
final, infinite excursion from the graph of Xf. The rest of the proof is devoted to that 
calculation. 
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Let U be the first time when the final, infinite excursion of B t from the graph of Xf 
hits the graph of Xf +s . Let 5\ = \X V — X v +S \. First, we will condition on Si. The process 
{B t ,t > U} is a Brownian motion conditioned not to hit the line t — > By + Si + fat. By 
subtracting the drift and flipping the process to the other side of the horizontal axis, we 
may consider a Brownian motion Bf starting from Si, with drift fa, conditioned not to 
hit 0. We will estimate the amount of the local time this process spends on the graph of 
a solution Y t to (1.1) with fa replaced by —{fa — fa), fa replaced by 0, and B t replaced 
by Bf. 

Let H Sl be the excursion law for excursions above the level Si for Brownian motion 
with drift fa, conditioned not to hit 0. Let Fx) denote the set of excursions with infinite 
lifetime. We will compute H Sl (Foo). Let Qp 2 be the distribution of Brownian motion with 
drift fa. Then 

H 5 i (Fqo) = lim - z • Q^+^Too < T Sl | < T Q ) 



lim 



1 Q% +z (Too < T Sl and < T ) 



,40 z g^(T oo <T ) 



= lim 



1 Q^ +z (Too<T 5l ) 



zioz g|+ 2 (T 00 <To)' 

Recall that the scale function S(y) for Brownian motion with drift fa is equal to 
exp(— 2fay/a 2 ). This gives 

H^)^ 1 S(h + z) - S(Sl) S M- S <°) 



zio z S(oo)-S{Si) S(6 1 + z)-S(0) 
= lim - e M-2fo(Si + z)/a 2 ) - exp(-2/3 2 5i/(7 2 ) - 1 



zio z - exp(-2/3 2 5i/a 2 ) exp(-2/3 2 (5i + z)/a 2 ) - 1 

= 2fa 

a 2 [l-exp(-2/3 2 5 1 /a 2 )]' 

If we fix arbitrarily small e > then for sufficiently small Si > we have 

^<^ 1 (^oo)<(l + £)^. (3.26) 
Si 6i 

We proceed to calculate the expected time to hit Si for Brownian motion with drift 
fa, starting from Si — z and conditioned not to hit 0, where z G (0, Si). If we take 



s(z) = exp 



- / 2fa/a 2 dy 
Jo 



exp(-2faz/a 2 ), 
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and 



S(z) = f s(y)dy = £-[1 -exp(-2/fcz/<7 2 )], 

then formula (9.9) on p. 264 of Karlin and Taylor (1981) yields 

E 5l -*(T 5l | T Sl < T ) 



2[S(6 1 )-S(6 1 -z)] 
S(5 1 )S(5 1 -z) 



f Sl ~ z S 2 (y) , f 1 

Jo v 2 s(y) Jst-z 



51 s^PM-sC*)] 



a 2 s(y)S(S 1 ) 



dy 



£[1- exp(-2/3 2 5 1 /a 2 )] - ^[1 - exp(-2/3 2 (5 1 - ^)/a 2 )] 



2/3 2 



^[1 - exp(-2/3 2 5 1 /a 2 )]^[l - exp(-2/3 2 (5 1 - ^)/a 2 ; 



x 



i g,[l-exp(-2/3 2?/ /a 2 )] 2 



+ 2 



/ a 2 exp(-2/3 22 //a 2 ) 

51 / ^[l-exp(-2/3 22/ /a 2 )] 



■ x 



Sl - Z ^ 2 exp(-2/3 2 y/a 2 )^[l -exp(-2/3 2 5 1 /a 2 )] 
^[1 -exp(-2^ 1 /a 2 )] - ^[1 -exp(-2/3 2?/ /a 2 )] 

^P2 Z P2 



2[exp(-2/3 2 (£ 1 - ^)/a 2 ) - exp(-2/3 2 ^/a 2 )] 
^[1 - exp(-2/3 2 5 1 /a 2 )][l - exp(-2/3 2 (5 1 - ^)/a 2 )] 



■ x 



Jo 



2\]2 



5 — i ^[l-exp(-2/3 2 y/ C r 2 )] 



+ 2 



/ 



a 2 exp(-2/3 2 y/a 2 ) 
*i ^[1 - exp(-2/3 2 y/a 2 )][exp(-2/3 2 y/a 2 ) - exp(-2/3 2 5 1 /a 2 )] 



dy. 



a 2 exp(-2/3 2? //(7 2 )[l - exp(-2/W 2 )] 
The expected lifetime of an excursion below Si for the Brownian motion with drift /? 2 , 
starting from Si and conditioned not to hit is therefore equal to 

lim-E <5l - 2 (T (5l \T 5l <T ) 

2 ^exp(-2/3 2 5 1 /a 2 ) f' 1 if^t 1 ~ exp(-2/3 2 y/a 2 )] 2 



z{0 Z 



202 



g-[l - exp(-2/3 2 5 1 /a 2 )][l - exp(-2/3 2 <f 



r 

i/a 2 )] 7o 



a 2 exp(-2/3 22 //(7 2 ) 



dy 
(3.27) 



for small <5i > and some C\ depending on /? 2 and cr 2 but not on 8\. 

Fix arbitrarily small e > 0. We are ready to derive estimates for the total amount 
of local time, say L^+, that spends on the graph of Y t . 

On one hand, the estimate (3.26) shows that is stochastically bounded by an 
exponential random variable with mean Si, for sufficiently small Si > 0. 
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Let Vg(a) be the time spent by £? t 9 between Y t and the horizontal axis before the 
time when L 9+ hits a. Since Y t is non-increasing, the estimate (3.27) can be used as an 
upper bound for the expected duration of an excursion below Y t , for every t > 0. Fix 
arbitrarily large b < oo and arbitrarily small e > 0. We have from excursion theory, 

Evjj"(6<5i) < WickJi, 
and so, for sufficiently small 8i > 0, 

POtf(Wi) > M < = — *i < e. 

die e 

We see that with probability greater than 1 — e, the distance between F f and the horizontal 
axis remains greater than Si — 8ie(fa — fa), at least until the time when L 9+ exceeds b8i. 
On this time interval and given this event, the intensity for the arrival process of the 
infinite excursion is bounded above by (1 +e)/(6i(l — e)), by (3.26). Hence, EL^ /Si can 
be made arbitrarily close to 1, by choosing large 6, then small e and finally small \8\ > 
(note that 8i<\8\). 

Finally, in order to obtain an unconditioned estimate for EL^+, we have to average 
over the possible values of 6±. Let v±(a) be the time when the local time of Bf (defined 
earlier in the proof of part (ii)) reaches a. The same argument which gives (3.22) yields 
the following estimate, 

EL 9 + = E8i + o(8) 

= \8\ -E(pi - p 2 )v A (a) + o(8) 

= \S\ - ^(1 - eM-MPi - fo)/* 2 )) + o(5) 

Pi — P2 

= \8\ exp(-2a(/?! - /3 2 )/a 2 )) + o(8). 
Adding this quantity to (3.25) gives the formula in Lemma 3.4 (ii). □ 

Lemma 3.5. Fix x,a, fii, fa > and assume that (3i — fa > 0. Then 

Var (L^ s | L^=a) = ^-^-(1 - exp(-2a(/3 1 - fa) /a 2 )) + o(8), 

Pi — P2 

for 8 I 0. The same formula holds if x < and 8 f 0. 

Proof. First suppose that x, 5 > and recall the notation and definitions from the proof 
of Lemma 3.4 (i). It follows from (3.23) that 

40 



We have using (3.22), 



E (< + (a) - a ? = E «(«) - ^^(a)) 2 (3.28) 
= E(S„ 8 4(a) ) 2 - 2ftE [B 8 V4{a) v 4 (a)] + E(v 4 (a)) 2 

= Ev 4 (a) - 2/3iE [sJ 4(a) U4(a)] + E(i> 4 (a)) 2 

(1 - exp(-2a(A - /3 2 )/a 2 )) + o(<S) - 2&E [a* ^(a)! + E(z; 4 (a)) 2 . 



Recall that v 4 (a) is bounded by 5/ (Pi — /? 2 )- Hence, 



E(v 4 (a)y<5'/(Pi-p 2 y, 



(3.29) 



and 



®B 8 Ma) v 4 (a) < (E(B 8 Ma) ) 2 E(v 4 (a)) 2 ) 

= (Ev 4 (a)E(v 4 (a)) 2 ) 1/2 
5 



1/2 



< 



Pi -Pi 
= o(d). 

Combining (3.28)-(3.30) yields 



(l-exp(-2a(A-/3 2 )/(7 2 ))+o(5) 



S 2 /(Pl~P2)'- 



(3.30) 



1/2 



E(L 4 + - a) 2 = -(1 - exp(-2a(A - /3 2 )/a 2 )) + 0(6). 



This implies 



ft- #2 



Var L 4 + , = E(L 4 + , - EL 4 + 

U4(aj \ 114(a) V4(a 



114(a) 1)4(0)' 



- a) 2 - (ELt+ a) - a) 



i) 4 (a) 



A -P: 



■(l-exp(-2a(/3 1 -/3 2 )/a 2 ))+o(5) 



" exp(-2a(/3 1 - p 2 )/a 2 )) + 0(6) 

Pi — P2 

■(l-exp(-2a(/3 1 -/3 2 )/a 2 ))+o(5). 



Pi-P: 



Similarly to the proof of Lemma 3.4 (i), we have Var (L^ s \ = a) = Var L 4 ^, which 
combined with the last formula proves the lemma in the case x > 0. 
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Now consider the case x, 5 < 0. We argue as in the proof of Lemma 3.4 (ii) that we 
have to add a contribution from the local time on Xf +s generated by the infinite excursion 
of B t below Xf . We have shown in the proof of Lemma 3.4 (ii) that the local time on Xf +S 
generated by the infinite excursion is stochastically bounded by an exponential random 
variable with mean 5 so its variance is bounded by 25 2 , and, therefore, the contribution to 
the variance from the infinite excursion is negligible. The same formula holds in the case 
x < as in the case x > 0. □ 



Recall that denotes the local time of B t — Xf at 0. The following result is 
analogous to Trotter's theorem on the joint continuity of local times for Brownian motion 
(see Karatzas and Shreve (1988) or Knight (1981)). 

Theorem 3.6. Assume that /?i,/?2 > and f3\ — f3 2 > 0. There exists a version of the 
process (x, t) — > which is jointly continuous in both variables. 

Proof. Note that Xf and Xf increase at the same rate when B t does not lie between Xf 
and Xf, and by the assumptions on fti and f3 2 , they grow closer together when B t does 
lie between them. Therefore, for all x,y and t > 0, 

\Xf-Xf\<\x-y\. (3.31) 

Define G(x) = EL^. The excursion law for Brownian motion below the line t — > f3 2 t 
gives mass 2f3 2 /& 2 to excursions with infinite lifetime, by Lemma 3.1 (iii). By excursion 
theory, the waiting time (in terms of local time) for the first excursion with infinite lifetime 
is exponential with mean a 2 /(2/3 2 )- This says that the distribution of is exponential 
with mean a 2 /(2/3 2 )- This and the strong Markov property applied at the first time when 
B t intersects Xf imply that for some c\ < oo and all x, we have G(x) < c\. An easy 
conditioning argument that combines this observation with Lemma 3.4 shows that for all 
x and y, 

\G(x) - G(y)\ < c 2 \x - y\. (3.32) 

The process (Xf , B t ) is strong Markov, and Lf is an additive functional. So by the 
Markov property, 

= G(Xf-B t ). 
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Therefore 

E[Z&-L? \r t ] < Cl . 

Also, using (3.31) and (3.32), 

|E[U& - LSo) - (Lf - L\) | ^]| = |G(Xf - B t ) - G(X» - B t )\ 

<c 2 \Xf-Xf\ 
< c 2 \x - y\. 

By Bass (1995), Proposition 1.6.14, 

E[sup|Lf -L y t \ 4 } < c 4 \x-y\ 2 . (3.33) 
t 

By Kolmogorov's criterion and standard arguments (cf. the proof of Proposition 1.6.16 of 
Bass (1995)), we deduce that there exists a version of Lf that is jointly continuous in x 
and t. □ 



A classical Ray-Knight theorem (see Knight (1981), Revuz and Yor (1991) or Yor 
(1997)) asserts, roughly speaking, that if Lf is the local time for the standard Brownian 
motion then x — > is a diffusion for certain stopping times T. As a part of that theorem, 
the infinitesimal parameters of the diffusion are also given. We prove a similar result for 
our family of local times, with T = oo. Recall that we assume that So = 0. 

Theorem 3.7. Suppose that p\,p2 > and (3\ — > 0. The distribution of 

is exponential with mean a 2 /(2/3 2 ). The process {L^ 7 x > 0} is a diffusion with the 

infinitesimal drift 

Ma) = -7r^Sr(l - exp(-2a(A - /3 2 )/a 2 )), 

Pi — P2 

and infinitesimal variance 

5 » = ir^ir^ - exp(-2a(A - /3 2 )/a 2 )). 
Pi — Pi 

The process {L^f, x > 0} is a diffusion with the infinitesimal drift 

Ha) = "TT^V + 7 T^Vexp(-2a(/3 1 -/3 2 )/a 2 ), 

Pi — P2 Pi — Pi 
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and the same infinitesimal variance 



^ 2 («) = ~n~~ird - exp(-2a(A - (3 2 )/a 2 )). 
Pi — Pi 



Proof. We have already shown in the proof of Theorem 3.6 that the distribution of 
is exponential with mean a 2 /(2/3 2 ). 

The Markovian character of the process {L^ 7 x > 0} at any fixed "time" x = y 
follows from the independence of the Poisson processes of excursions of B t below and 
above X\. The same remark applies to {L^f, x > 0}. The infinitesimal parameters of the 
processes were calculated in Lemmas 3.2 and 3.3. 

The process x — > is continuous, by Theorem 3.6. Since its infinitesimal drift is 
bounded and the infinitesimal variance is nondegenerate, there is a unique (in law) Markov 
process with this infinitesimal drift and variance (cf. Bass (1997), Section IV. 3), and this 
Markov process is in fact a strong Markov process. □ 

Theorem 3.8. Suppose fli, f3 2 > and /?i — f3 2 > 0. For fixed t > 0, we have a.s., for all 

x, xi, x 2 e ffi, 

= exp(-2Ln/3i-/3 2 )/a 2 ), 



d 

dy 

and 



-X? 



X^ - X? = / exp(-2Ln/3i - (3 2 )/a 2 )dx. 

J Xl 



Proof. First we will prove an estimate analogous to (3.22) except that it will hold for 
v^(a) itself rather than its expectation. Recall the notation and definitions from the proof 
of Lemma 3.4 (i). 

The following estimate is completely analogous to (3.17) except that we state it for 
the process Bj rather than Bf : so Si is replaced by 5 in the bound. 

Cn<J 



¥ y (v 7 (s (l-e)) > u ) < 



e 2 so<j 2 



We can further modify the estimate by replacing sq(1 — e) with so, so that 

cn<y(l - e) 



^(v 7 (so)>u /(l-e))< ,. , 

£ oqO 
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This and (3.16) imply that for small e and 5 we have 



v 7 {s ) < u /(l -e)< 2s 5(l + e) 2 /a 2 (3.34) 

with probability greater than or equal to 1 — C\\8{1 — e)/(e 2 so<j 2 ). The inequality (3.34) 
is analogous to (3.11) and can be used in the same way as in the argument between (3.19) 
and (3.20) to prove a formula analogous to (3.20): 

vi(a)=v 4 (js ) < -(l-exp(-2a(/3i - f3 2 )(l + ef /a 2 )), (3.35) 

Pi — P2 

where 5 in i>|(a) indicates the dependence of vf(a) on 5. The above argument requires 
that we can use an estimate analogous to (3.34) at every stage of the inductive procedure, 
i.e., at every stopping time v^imso) for m = 1, 2, . . . ,j — 1. All of these estimates hold 
simultaneously with probability greater than 

l-(j-l)c 11 S(l-e)/(e 2 s a 2 ). 

This shows that the probability that (3.35) fails to hold is smaller than 

(j-l)c 11 5(l-e)/(e 2 s a 2 ). 

Now fix arbitrarily small e± > and let 8 k = (1 — Si) k . Let A k be the event in (3.35) with 
5 replaced by 8 k , i.e., 

A k = |^ fe (a) < ^-^-(l-exp(-2a(/3 1 -/3 2 )(l + £ ) 3 /a 2 ))}. 

We have 

oo oo 

- l)cn4(l - e)/(e 2 s a 2 ) = ]T(j - l)c n (l - e 1 ) k (l - e)/(e 2 s a 2 ) < oo, 

k=0 fc=0 

so only a finite number of events A k may fail to hold. Consider an u and ko such that 
all events A k , k > k 0l hold for this u>. Suppose that 8 £ (0, 8k ). Then 5 G [S^-i, 8 kl ] for 
some ki > ko. Since A kl holds, we have 

v S 4 (a) < v^(a) < ^V(l -exp(-2o(A - +ef/a 2 )) 

Pi — P2 

< 5 f (1 ~ exp(-2a(A - + ef/a 2 )). 
Pi — Pi 
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This inequality holds with probability one for all sufficiently small 5 > 0. Since e > and 
E\ > are arbitrarily small, we see that a.s., 

limsup < TT^-Tra - eM-MPi ~ fc)/* 2 )). 

The same lower bound can be obtained for liminf in a completely analogous way, so with 
probability one, 

lim 4(a)_ = l-exp(-2a(/3 1 -/3 2 )/(7 2 ) 
5^0+ 5 p\ - (3 2 

Suppose a > and let v(a) = inf{t > : = a}. Fix some x G K. and consider 

S > 0. We will first find the right hand side derivative jj-X^,^ . Let T = inf{t : S t = 

y \ I y =X 

Xf } and let £/i be the amount of time spent by B t between the graphs of Xf and X^ +s 
on the time interval [0, T]. We will write U 2 to denote the amount of time spent by Bt 
between the graphs of Xf and X^ +s , between times T and v (a). 

If x > then t/i = 0. If x < then Z7"i is not greater than the amount of time Us 
spent by B t between the lines t — > x + and t — > x + 5 + Pit, until the hitting time T. 
Standard arguments show that for any arbitrarily small e > 0, we have Us/5 2 ~ e — > as 
5^0, a.s. Note that the distance between X* +s and decreases by (f3\ — p2)« on any 
interval where the Brownian motion B t spends u units between these functions. Hence, 

Kta) - X v(a) =8-Vh- + U 2 ). 

The random variable U 2 may be identified with vf(a), so (3.36) gives for any fixed a, a.s., 



dy v ^ 



yx+S -yx 

,. via) ^v(a) 

= lim — 

y=x 

= ]im S-(Pl-lh)(U 1 + U2) 



5^0 5 

= 1 - lim (fl-fo)^ _ lim (Pi-^Hia) 

= 1 - - (1 - exp(-2a(A - (5 2 )/a 2 )) 
= exp(-2L- (a) (A-/3 2 )/a 2 ). 

The above holds simultaneously for all rational a, with probability one. Since t — > and 
t — > — XI are continuous monotone functions, an elementary argument can be used to 
extend the last formula to fixed times, i.e., 

d+ 



dy~ Xt 



= exp(-2Ln/3i-/3 2 )/^), (3.37) 

y=x 
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simultaneously for alH > 0, a.s. 

Fix some t > 0. By Fubini's theorem, (3.37) holds for almost all x, a.s. We have 
|Xf — X t 2 | < \y — z\ for all y and z. Since the function y — > Xf is Lipschitz, it has a 
derivative almost everywhere and so for a fixed £, we may replace the right hand derivative 
with the usual derivative in (3.37), for almost all x. The function x — > is continuous, so 
the derivative in (3.37) is equal almost everywhere to a continuous function. This implies 
that the derivative is equal to the function everywhere. This proves the first assertion of 
the theorem. The second one follows from the first one and from the Lipschitz character 



Remark 3.9. Suppose that Xf are solutions to (1.1) and assume that Pi, 02 > and 
Pi — Pi > 0. Fix some t > and consider the function y — > Xf . We will sketch an 
argument showing that y — > Xf is C 1+1 for every 7 < 1/2, i.e., that the function has a 
derivative which is Holder continuous with Holder exponent 7. 

Fix any zGt. With probability 1, Bt 7^ X*, and with strictly positive probability, 
there exists e > such that B s 7^ X^ for all y G (-2 — e, z + e) and s > t. It follows that 
if a local property holds for the function y — > with probability 1, it must hold for 
y — > L^, with probability 1. Since y — > is a diffusion, its paths are Holder continuous 
with exponent 7 for every 7 < 1/2. It follows that the same is true of y — > Lf. Theorem 
3.8 now implies that y — > is C 1+7 for every 7 < 1/2. The same argument shows that 
y -> Xf is not C 3 / 2 . 

4. Time and direction of bifurcation. We will first address the question of the direction 
of bifurcation for the equation (1.3). We will say that a positive bifurcation occurs if for 
some ti we have X t > B t for all t > t±. The definition of a negative bifurcation is 
analogous. If f3\ and #2 have the same sign then it is easy to see that a bifurcation will 
occur with probability one and its direction will be the same as the sign of /Vs. If /?i > 
and < then there will be no bifurcation. The next theorem deals with the only 
remaining, non-trivial case. 

Theorem 4.1. Consider the equation (1.3) with t = x = 0. Assume that Pi < and 
p 2 > 0. Let 

j ~ a 2 /K+i)r(l/(a J + 1)) ' 
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for j = 1,2. The probability of a negative bifurcation is equal to Ai/(Ai + A 2 ). When 
«i = «2 = 0, tie formula simplifies to \ + | /?2 1 ) - 

Before proving Theorem 4.1 we present a lemma which may have some interest of 
its own. 

Lemma 4.2. Assume that Pi < and 02 > 0. Consider a solution Xt to (1.3) with 
to = x o = 0. There exists 7 > 0, depending on a±, «2, such that Xt/t 1 ' 2+1 converges 

in probability to as t — > 0. 

Proof. Let us assume that — 1 < cti < < a^- The other cases may be treated in a similar 
way. Let U = sup{s < t : X s - B s > -t 1 ' 2 }. For s e (U, t) we have X s - B s < -t 1 / 2 so 
for such s, |dX a /ds| < /5it (1 / 2 > Ql . It follows that 



*t - ^/| < -(t - t/)/3it (1/2)Ql < -/3it 1+(1/2)Ql , 



and so 



X t < (B t - B v ) + (B v - X v ) + (X v - X t ) 




This implies that 



E\B t -X t \ ai l 



{B t -X t >0} 




< Clt (l/2)ai +C2 | /3l |a lt a 1 + (l/2)a?_ 



Recall from (2.2) that 








From this we have the following estimate 
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Since ct\ > — 1, the exponents 1 + {l/2)a\ and 1 + ai + (l/2)ct 2 are greater than 1/2 and 
so for some 7 > and every C5 > 0, liminf^o c^EXt/t 1 / 2 ^ > 0. It follows that 

limP(X t A 1 /2+7 < _ C6)=05 

for every cq > 0. 

Recall that a 2 >0. Since 

x t < I [3 2 \x s -B s rn {Xs _ Bs>0} ds, 

Jo 

an elementary argument shows that for small t we have X t < 2f3 2 t if B s < 1 for all 
s G (0, t). It is clear that P(max s€ ( t ) B s > 1) goes to as t — > so 

\im¥(X t /t > 2f3 2 ) = 0. 

This and (4.1) prove the lemma. □ 



Proof of Theorem 4.1. The assertion of the theorem deals only with probabilities, so 
we can use any solution to (1.3), as we have uniqueness in law by Theorem 2.1. The same 
theorem shows that a solution X t may be constructed so that (X t , B t ) is a strong Markov 
process, and hence we may apply excursion theory to it. Recall the discussion at the 
beginning of Section 3. The same analysis of excursion laws and the exit system applies 
to the solutions of (1.3) for arbitrary a\,a 2 > —1. Let us briefly recall the facts that we 
will need in our present argument. Let D = {(6, x) G M 2 : b = x} and let (iiP, dL) be an 
exit system for the process of excursions of (B u X t ) from the set D. The generic excursion 
may be denoted (ej, e 2 ). By the translation invariance of the Brownian motion B t and the 
equation (1.3), the distribution of {e\ — x, e 2 — x) under H x is the same for every iet 
Let this distribution be called H\. Let H\ + denote the part of the measure Hi which is 
supported on excursions with e\ > e 2 and let be the part supported on the set where 
e] < e 2 . Let H 2+ be the distribution of {e\ — e 2 , t G (0, z/)} under Hi + and let H 2 - have 
the same definition relative to i!Zi_. We have, up to a multiplicative constant, 

H 2+ (A) = lim ^Ql(A), (4.2) 

xj.0 \x\ 

where Q\ stands for the distribution of the diffusion Y t with the same infinitesimal variance 
as Brownian motion (i.e., a 2 ) but with drift —/3i\Y t \ ai : killed at the hitting time of 0. We 
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will normalize H 2+ as in (4.2). We will have to prove that the following formula gives the 
correct normalization for H 2 _, 

H 2 _(A) = \im-^-Q x _(A). (4.3) 

Here Q x _ denotes the distribution of diffusion Z t with Brownian quadratic variation (name- 
ly, a 2 ) and drift — /3 2 |Zt| a2 , killed at the hitting time of 0. 

The proof that (4.3) is the correct normalization for H 2 - can proceed exactly as 
the proof of Lemma 3.1 (iv), thanks to Lemma 4.2. It only remains to find and compare 
the formulae analogous to those for H 2+ (F h ) and H 2 -(F h ). Recall that F h is the event 
that the difference between the maximum and the minimum of an excursion exceeds h. 
The scale function for a diffusion on (0, 00) with infinitesimal drift /j(x) = —f3ix ai and 
variance a 2 is given by (see Karlin and Taylor (1981), p. 194), 



S(x) = I exp ^ — J ~^' 2 ~ / dzj dy = J exp ( - / — dz ) dy 



By (4.2), 



H 2+ (F h ) = lim-Ql (T h < To) = lim - ■ S J^ §g 
x.10 x + xio x S(h) - S(0) 

, 1 £«p(§8g=TO)* 1 

= lim — • 7 ( = 7 r . 

If we use (4.3), we obtain in the same way 

H 2 -(F h ) = 



which implies that 



lim H 2+ (F h )/H 2 _(F h ) 



and this confirms that the normalization in (4.3) is correct. 

The probability for the process Y t starting from 8 never to hit is equal to 



lim 

6— >oo 



s(S) - s(o) _ So fq^ ggrgO j; 
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It follows that H 2 +(F 00 ) : i.e., the measure given to positive excursions which do not return 
to is given by 



lim- 

5^0 5 



J °exp thrill 



<r 2 K+i) 



dy 



/o"«p(i^)* 



exp 



ai + l 



a 2 (ai + l) 



-i -i 



( 



-Mi 

cr 2 (a! + l) 



l/(ai+l) 



r(V(«i + i)) 



(ai + 1) (-2ft)i/(«i+i)( ai + df 

- — Ai. (4.5; 



a 2 /(«i+i)r(l/(a 1 + 1)) 



df 



An analogous formula holds for A2 = H 2 -(F 00 ). The processes of excursions on both sides 
are independent so the probability of the negative bifurcation is the same as the probability 
that the first arrival of an infinite excursion in the Poisson process on the negative side 
comes before the analogous event on the other side. The probability in question is the 
ratio of Ai and Ai + A2. □ 



Remark 4.3. Mike Harrison pointed out to us that Theorem 4.1 may be proved without 
using excursion theory. One can calculate the probability that the diffusion X t — B t will go 
to infinity using an explicit formula for the scale function of this diffusion. The excursion 
theory approach has its advantages, though. First, excursion theory seems to be the right 
tool for the proof of Theorem 4.4 below. Second, the excursion theory may be used to 
find the positive bifurcation probability when the vector process (X t , B t ) is Markov but 
X t — B t is not. The solution of (1.2), studied in Burdzy, Frankel and Pauzner (1998), is 
an example of such a situation. 

Let T* denote the bifurcation time, i.e., let T* be the supremum of t with X t = B t . 

Theorem 4.4. Consider the solution to (1.1) with to = xq = 0, (3\ < and 02 > 0. Then 



Proof. By Remark 3.2 (ii), the distribution of the excursion law on excursions with finite 
lifetime remains the same if we change f3 to — (3. Hence, the formula (3.1) applies in the 

case Pi < and p 2 > 0, and we have 
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for V s — inf{t > : L t > s}. By excursion theory, the infinite excursion of B t from X t ° 
occurs independently from finite excursions in the Poisson point process of excursions, so 
the expected bifurcation time is equal to 



oo 

AST 



ET* = / Xe~ MV s ds, 



where A is the intensity of the Poisson process arrival for infinite excursions. We have 

2(|/3i| + |/3 2 |) 



A = 



a 2 



from (4.5), taking into account infinite excursions on both sides. It follows that 

f°° r°° 1/11 

ET* = / \e~ Xs EV s ds = / Ae" As — + — ) sds = - — + 



o Jo 



\0i\ \fo\J A VIA | \h 



^ < 1 1 \ a 2 

+ 



2(|/5i| + 1^|) \fo\J \2W 



Remark 4.5. (i) It is also the case that 



ET, 



if > 0. We leave the proof to the reader. 

(ii) A similar result can be obtained for any values of «i, a.2 > — 1 but the formula 
does not seem to have a compact form, so we only sketch how it can be obtained. The 
proof of Theorem 4.4 needs two ingredients. One of them is the expected amount of local 
time before the infinite excursion occurs. This is equal to the expectation of the minimum 
of two independent exponential random variables whose expected values are inverses of the 
quantity in (4.5) (for (ai,Pi) and (a^,/^))- 

The second ingredient is the expectation of the inverse local time at s, for the 
process with finite excursions only. This is equal to s times the expected lifetime of a 
finite excursion under the excursion law. Here is how we can calculate this quantity. For 
arbitrary ol\ > —1 we write as in (4.4), 

s(x) = exp 



a 2 (ai + 1) 
f 2/3iy ai+1 



S(x)= I exp(- g2(Qi + i) )dy. 
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For < x < y < oo, the Green function for Brownian motion Y t with drift —j3\ (the 
negative sign is due to restriction of the excursion law to finite excursions) is given by (see 
Remark 3.3 on p. 198 of Karlin and Taylor (1981)), 

2[S(x)-S(0)][S(oo)-S(y)] 
a*s(y)[S(oo) -S(0)} ■ 

Hence, the expected lifetime of an excursion is equal to 

WxJ 0{x ' y > dy l Z$, x ] <rM»)[S(°°) - S(0)] S 



» ! [S(»|-S(0|]i, »(y) 



Jo s (y) 



«1+1 



2/3iz 
cr 2 (ai + l) 



dz 



a j Q exp ^ o-2( ai+1 ) J az v exp ^ 0-2(^+1) y 

Adding this to the analogous quantity for oii gives the expected lifetime of a finite excursion 
under the excursion law. 

5. Lipschitz approximations. In this section we will address the question of how well 
a Lipschitz function can approximate a Brownian path. Our analysis will be based on the 
fact, proved in Lemma 5.2 below, that a certain solution to (1.1) may be looked upon 
as a Lipschitz approximation to B t . 

Our first lemma consists of two elementary observations which are designed to help 
develop the mental picture of the solutions Xf of (1.1), in preparation for Lemma 5.2. 

Lemma 5.1 Let Xf denote the solution of (1.1) with Xq = x. 

(i) Ifx<y then Xf < Xf for all t G R, a.s. 

(ii) For a fixed t, the function x — > Xf is continuous a.s. 

Proof, (i) Suppose that we have x < y and Xf — X^ for some sGl. The two functions 
Xf and X\ are not identical since Xq = x ^ y = X%, but they are both solutions to (1.1) 
with t = s and xq = Xf. This contradicts the uniqueness of solutions to (1.1). 

(ii) Consider any sequence x n converging monotonically to x^ G K. By (i), the 
sequence Xf n is also monotone in n, and by the Lipschitz property it must converge to a 
limit Xf°. The Lipschitz property of the Xf n, s implies that of Xf° . We can show that 
Xf° is a solution to (1.1) using the same argument as in the proof of existence in Theorem 
2.2 for (1.1). We must have X^° = Xq°°, so the uniqueness of the solutions to (1.1) implies 
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that Xf° = Xf°° for all t, a.s. We have shown that x n -> implies Xf n -> Xf°°. This 
completes the proof. □ 

Lemma 5.2. Assume that fli < < fa. For almost every oo there exists a unique x = x(u) 
such that the solution Xfi to (1.1), that is, the solution satisfying Xfi (uj) = x(uj), has the 
property that there exist arbitrarily large t with Xf(uj) = B t {u)). 

It is easy to see that if fa < < fa then with probability 1, all solutions Xf have 
the property that there exist arbitrarily small t > — oo such that Xf = B t . Lemma 2.13 
shows that a result analogous to Lemma 5.2 holds when fa < < fa, and we require that 
the solution intersects the Brownian path for arbitrarily small t > — oo. 

Proof. We will first prove the existence. The law of the iterated logarithm easily implies 
that for some random x > 0, the functions t — > x + fat and t — > — x + fat stay above and 
below the trajectory of B t , for t > 0, resp. This shows that there exist both large and 
small (random) x such that Xf does not intersect the trajectory of B t for t > 0. 

Let A be the set of all x such that Xf > B t for all t greater than some ti = t\{x). 
By Lemma 5.1 (i) and the above remarks, the set A is a non-empty semi-infinite interval. 
We will show that it is open. Consider an x such that Xf > Bt for all t greater than some 
£i. Then Xf = Xf x + fa(t - h) for t > t x . Let c x = mf{Xf - B t : t > h + 1} and note 
that ci > 0, by the continuity of Xf — Bt. By Lemma 5.1 (ii), the function y — > Xf i+1 is 
continuous so we can find e > such that Xf i+1 > Xf i+1 —c\/2 for all y > x — e. It follows 
easily that for such y, we have X\ = Xf i+1 + fa{t — 1\ — 1) and so Xf > B t for t > t\ + 1. 
This proves that A is open. The same is true of the set A' of x's with the property that 
Xf < B t for all t greater than some t\ = t\(x). Hence, (A U A') c is non-empty and so we 
must have at least one x for which Xf = B t for arbitrarily large t. 

We turn to the proof of uniqueness. Suppose that with positive probability there 
exist x\ < X2, such that both trajectories Xf 1 and Xf 2 intersect B t for arbitrarily large 
times t. Then we can find 6 > and p > such that with probability greater than p, there 
exist xi and X2 with xi > x\ + 5 and such that Xf 1 and Xf 2 intersect B t for arbitrarily 
large times t. We will show that this assumption leads to a contradiction. 

Fix some 7 G (1/2, 1). By the law of the iterated logarithm, we can find a large ti 
with the following property. For every x, if \Xf \ > tj then X} 7^ B t for all t > t\, with 
probability greater than 1 —p/8,. 
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Consider solutions Xf and Xf to (1.1) with Xg = yi = t\ and X y t 2 = y 2 = -tj. 
We enlarge t±, if necessary, so that the event {|_B fl / 2 | > tj} U {\B tl \ > tj} has a probability 
smaller than p/8. If the event {|_B tl / 2 | > tj} U {\B tl \ > tj} does not occur and t± 
is sufficiently large, then both processes Xf 1 and X\ 2 must intersect the trajectory of 
B t between t±/2 and t\. Let T = sup{t < t\ : X\ x = B t }. Then the process {Y t = 
B T _ t — X^_ v t > 0} is a Brownian motion with drift f3 2 if Y t < and /?i if F t > 0. Note 
that the distribution of the process Y t does not depend on t\. We can apply Proposition 3.3 
to the local time Lj of Y t at 0, to see that Lj /t — > (l/|/3i| + 1 / 1 /?2 1 ) 1 = A, as t — > oo, a.s. 
Enlarge ti again, if needed, so that > \t\/A with probability exceeding 1— p/8. Let ?y 
be the expected number of positive excursions of Y whose height does not exceed 6, whose 
duration exceeds 1, and which start at a time t with Lj < 1. Then the total number of such 
excursions which start at times t with Lj < Xti/4 has a Poisson distribution with mean 
rjXti/A. We make t\ large enough so that with probability greater than 1 — p/8, the total 
number of such excursions which start at times t with Lj < Ati/4, is greater than r]Xti/8. 
Collecting the above facts, we see that with probability greater than 1 — 3p/8, we have all of 
the following: (i) T G (ti/2, £i), (ii) the local time for the process Xf 1 —B t at accumulated 
between times and t\ exceeds Ati/4, and (iii) the number of negative excursions of 
Xf 1 — Bt whose absolute height is less than 5, the duration is greater than 1, and which 
lie within interval (0, is greater than rjXti/8. Note that dXf 2 /dt — dXf 1 /dt = f3 2 — Pi 
for any t within such an excursion provided Xf 1 > X\ 2 + 5. The last observation shows 
that if X,f > X% 2 + 5 then Xf^ > Xf 2 + 5 + 2(f3 2 - p 1 )r]Xb 1 /8. If h is sufficiently large 
the last inequality cannot hold because we would have tj > —tj + 5 + 2(f3 2 — /3i)rjXti/8. 
We conclude that with probability greater than 1 — 3p/8, we have X^ 1 — Xq 2 < 5. 

We reformulate the last statement in terms of x\ and x 2 . Using Lemma 5.1 (i), 
we see that the probability that there exist x\ and x 2 with x 2 > x\ + 5, \X^\ < tj 
and \Xt 2 \ < tj is less than 3p/8. An earlier argument showed that the probability that 
x 2 > x± + 5 and \X^\ > tj or \X* 2 \ > tj is less than 2p/8. We conclude that the 
probability of x 2 > x\ + 5 is bounded by 5p/8, which contradicts our assumption. □ 

Consider equation (1.1) with —j3\ = (3 2 = (3 > 0. Let X£ denote the solution of 
(1.1) constructed in Lemma 5.2. That is, X? = Xf . 

Lemma 5.3. We have with probability 1, 

X* - B t X* - B t B t -X* Bt-Xl a 2 
hm sup = lim sup = lim sup — = urn sup > — . 

t^-oo log* t^oo log* t^-oo log* t^oo log* 2/3 
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Note the lim inf's as t — > oo are zero as X£ crosses B t for arbitrarily large t. 

Proof. Let X t be a solution to (1.1) with —f3\ = fa = —f3 and let Y t = B t — X t . 
The process Y t is a diffusion which spends zero time on the real axis, which behaves like 
Brownian motion with drift (3 when Yj < 0, and it is a Brownian motion with drift — (3 
when Y t > 0. By Karlin and Taylor (1981), Chapter 15.5, (5.34), the process Y t has a 
stationary probability distribution with a density 

w{y) = -o exp 



o ^"-f 1 9 

u \ u 

Let G M} be the process which has density ip(y) for every fixed t, and which has the 

transition probabilities of Y t . Let 

X t = [ sgn(Y s )(3ds 
Jo 

and 

B t = Y t -Y + X t . 

It is easy to check that B t is a Brownian motion with B = 0, and that X t solves (1.1) 
with B t replaced by B t and (3\ = —fii = f3. Moreover, X t has the property that X t = B t 
for infinitely many arbitrarily large negative and arbitrarily large positive t. If we now 
time- reverse B t and X u we will obtain a Brownian motion and a corresponding solution 
to (1.1) which satisfies the defining properties of X% . Hence, we may construct B t and the 
corresponding process X? by letting B t = B- t and X£ = X- t . 

The scale function S(y) for Brownian motion with drift —f3 is given by S(y) = 
exp(2[3y/a 2 ) (Karlin and Taylor (1981) Chapter 15.4). Let T a be the hitting time of a by 
the process Y. The mass H(Fh) given by the excursion law for the process Y t to positive 
excursions whose height exceeds h is equal to 

, m \ , 1 5(e) -5(0) , 1 exp(2(3e/a 2 ) - 1 
lim -P (T/i < T = lim - • -^4 = lim F ^ o ; 

eio e y U) eio e S{h) - 5(0) e[o e exp(2(3h/a 2 ) - 1 

- M 1 

~ ' exp(2[3h/a 2 ) - 1' 

Fix some small e > and let /ifc = /clog 2 • (1 — £)a 2 /(2,9). Let L t denote the local 
time of Y t at with L = 0, and let Afc denote the event that there exists a positive 
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excursion of Y t whose height exceeds h^, and which starts at a time t such that 2 fc < L t < 
2 k+1 . The probability of Ak is the probability that a Poisson random variable with mean 
Afc = 2 k H{F hk ) takes a non-zero value. Thus, P(A£) = e~ Afc . For large k, 

= 2* . ^ exp C-24 ■ HOS2 'P-^ 2 ) = 2» 4 ■ 2-*C-> = 2? • 2-. 

This implies that J2k ^O^fc) = Sfc e_Afe < 00 • By the Borel-Cantelli Lemma, only a finite 
number of the events A c k occur. Hence, 

limsup — ^— > limsup sup i — ^— : L t £ [2 fc , 2 fc+1 ] 
t^oo logL t fc^oo [logL t 

> limsup sup I — ^xr : A G [2 fc , 2 fc+1 ] 

fc^oo [log2 fc + 1 

. v ^ 

> limsup 



fc^oo (fc + 1) log 2 

Hog 2 • (1 -eW 2 

= hm sup — — — 

k ^J 2/3(A; + l)log2 

2/3 ' 

Since £ > is arbitrarily small and, by Proposition 3.3, lim^oo L t /t = (3/2, a.s., we obtain, 
with probability 1, 



Yt a 2 



lim sup 
t^oo log* 2/3 



A similar argument yields, 



— lim inf — — = lim sup — — = — lim inf — — = — . 

t^oo log* t^-oo log* t^oo log* 2p 

Recall from the first part of the proof that Y t = B_ t — X*_ t — Y . This combined with the 
results for Y t implies the proposition. □ 



The function t — > a + (3\t\ is Lipschitz with constant (3. For some random a, this 
function is greater than B t for every t, by the law of the iterated logarithm. Since the 
infimum of an arbitrary family of Lipschitz functions with constant (3 is again a Lipschitz 
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function with constant /?, there exists a smallest Lipschitz function with constant (3 
with the property that Z^ > B t for all t. Let Z~[ be the largest Lipschitz function with 
constant f3 such that Z~[ < B t for all t. Note that Z t + and Z^~ are not measurable with 
respect to a{B s , s <t}. 

Lemma 5.4. Assume that B = 0. We have with probability 1, 

hm sup — = hm sup — < — . 

t^-oo log* t^oo log* 2/3 



Proof. Consider ai,a2 > 0. Let 

A ++ = {3t > : B t = ai + /%}, A + _ = {3t > : S t = -a 2 - 
A_ + = {3t < : B t = ai - /%}, A__ = {3t < : B t = -a 2 + fit}. 

The probability that S t ever hits the line t — > a\+(3t is equal to exp(— 2aif3/a 2 ) (Karlin and 
Taylor (1975), p. 362). The probability that B t crosses the line a\ + fit at some t\ > and 
then crosses the line —a 2 —/3t for some t > t\ is bounded by exp(— 2aif3/a 2 ) exp(— 2a 2 /3/a 2 ), 
by the strong Markov property applied at £i. The probability of crossing first — a 2 — f3t 
and then a± + j3t is bounded by the same quantity. Hence, 

P(A ++ n A+_) < 2exp(-2(ai + a 2 )(3/a 2 ). 

The same estimate holds for ¥(A fl A_ + ), by symmetry. We obtain 

F(A ++ n A__) = P(A_+ n A+_) = exp(-2(ai + a 2 )/3/a 2 ), 

from the independence of the processes {B t ,t > 0} and {B t , t < 0}. It follows that 
¥(Z+ -Bo> ai ,B - Z~ > a 2 ) = ¥(Z+ > a u Z~ < -a 2 ) 

< F{[A ++ n A+-] u [A— n u n u n A+-]) 

< 8exp(-2(ai + a 2 )/3/a 2 ). 

Choose e G (0, 1). Let m > 8 be an integer large enough that (m— l)/(m(l— e)) > 1. 
We have for any y > 0, 

P(Z+ - Z " > y) < ¥ ( Z o ~ B o> jy/m, B - Z^ > (m - j - l)y/m) 

3=0 

( -2{m-l)y(3 \ 

< 9m exp - . 

\ ma 2 J 
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Fix some large b < oo. Consider an integer k > 0. Let n be the integer part of 



2/3(1 -e)b2 k 
cr 2 Hog2 



and let Xk = b2 k /n, and t h , = j2 k /n. We have, 



P(Z+ - > x k ) = ¥(Z+ - Zq > x k ) 



t k t 

3 3 



\ mo £ J 

( (m-l)2f3b2 k \ 

— 9m exp - 

\ mna 2 J 

( (m-l)2f3b2 k a 2 k\og2 

- 9m6X H m2/?(l- e )62*<r» 

= 9m exp ( -fc(m -l)log2 \ = ^ 2 _ fc(m _ 1)/(m(1 _ e „ 



For some ci < oo, using (m — l)/(m(l — e)) > 1, we obtain, 



E E p(^t-^>^)<ciE 2 ^- 9w - 2 " fc(m " 1)/(m(1_e)) 

fc=l 0<t'?<2 fc J ^ fe=l 

< Ci y 2/3( 9 1 ~ £)62fc • 36m • 2 -*<™- 1 >/<™< 1 -*» < oo. 
(j 2 /clog2 



fc=i 



By the Borel-Cantelli lemma, for all sufficiently large and all t k G [0, 2 fc ], we have 
Zj - ^ < x fc . If Z+ - Z~ < x k then for t G [(tj + *J_ x )/2, (*J + t} +1 )/2], 

3 3 3 3 

Zt - Zt < x k + 2(3\t-tj\ < b2 k /n + f32 k /n = (b2 k /n)(l + 0/b) = x k (l + 0/b). 

This implies that for large k, we have for all t G [0, 2 fc ], 

Z+-Z t " <x k (l + P/b). 
59 



We obtain 



z~y — z+ Z^~ — z+ 

hm sup — < lim sup sup — 

t^oo lOgl fc^oo t e [ 2 fe-i,2 fe ] l°g* 

zf - Zt 

< lim sup sup 



Xjfc(l + /?/&) 



< lim sup 

fc^oo 

< lim sup 



fc^oo (A; -1) log 2 

(b2 k /n)(l + (3/b) 
fc^oo (fc - 1) log 2 

62 fc (7 2 Hog2(l + (3/b) 



" ™T 2/3(1 -e)62*(fc-l)lo g 2 

_ ^ 2 (l + /3/^) 
2/3(1-6) ' 

Since e may be chosen arbitrarily small and b may be chosen arbitrarily large, with prob- 
ability 1, 

Z+-Z t ~a 2 

hm sup — < — . 

t^oo log* 2(3 

The result for * — > — oo follows by symmetry. □ 



Theorem 5.5. (i) With probability 1, 

hmsup — = hmsup — 

t^-oo log* t^oo log* 

= hm sup — : = hm sup — 

t^-oo log* t^oo lOgt 

Zt-Zt r Zt-Zt a 2 

= hmsup = hmsup — = — . 

t^-oo log* t^oo log* 2(3 

(ii) E\B t -X*\ = \- a 2 / (3, for every t G R. 

(Hi) E(Z+ - B t ) = E(S t - Zf) = | • a 2 / (3, for every t G R. 

Theorem 5.5 (i) shows, in a sense, that Z + and Z~ are as good Lipschitz approxi- 
mations to B t as X*. However, the comparison comes out differently when we look at the 
averages presented in (ii) and (hi). 

Proof. It is elementary to check that we always have Zf < X£ < Zf . This and Lemmas 
5.3 and 5.4 yield (i). 
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Recall the stationary density ip(y) for Y t from the proof of Lemma 5.3. This is the 
same as the density for the distribution of B t — X£. Hence 

m - X t | = \y\^ exp (--5- )dy=--j, 

which proves (ii). 

For a > 0, the probability that B t crosses the line a + f3t for some t > is equal to 
exp(-2af3/a 2 ) (Karlin and Taylor (1975) p. 362). This is the same as the probability of 
crossing the line a — fit for some t < 0. The probability that none of these events happen 
is [1 - exp(-2a/3/a 2 )] 2 , and so 

F(Z+ < a) = [1 - exp(-2a/3/a 2 )] 2 . 

This yields 



We similarly have EZ ( 



¥7+ - 3(j2 

= — 3o" 2 /(4,9), and by translation invariance, for every t, 
E(Z+-B t ) = E(B t -Zr) = ^. □ 



If we let a,\ = «2 = 1 and choose suitable (3i and in (1-3), then Y t = X t — B t is 
an Ornstein-Uhlenbeck process. Results for such a process, closely related to Theorem 5.5 
(i), can be found in the paper of Darling and Erdos (1956). 

Corollary 5.6. For any random Lipschitz function g(t) with constant f3 we have with 
probability one 

limsup ^"^ >a 2 /(4/j). 
t^J log* " /Kf) 



Proof. Suppose that X£ — B t = a for some t and a > 0. Let s be the largest time less 
than t such that B s = X*. Then we see that the quantity sup s<u<t \g(u) — B u \ cannot be 
smaller than a/2 for any Lipschitz function g(u) with constant (3, by comparing g(u) with 
the function u — > B s + a/2 + (u - s)j3. Since limsup^^ (X t * - B t )/\ogt = cr 2 /(2/3), for 
any Lipschitz function g(t) with constant /3 we must have 

limsup ^?"^ >a 2 /(4/3). □ 
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Corollary 5.6 sheds some new light on an old problem about strong approximations. 
Let us assume that a 2 = 1, i.e., we will consider now only standard Brownian motion. 
Suppose that {Vk}k>i are i.i.d. random variables such that \Vk\ < 13, a.s. Let S n = 
Ylk=i Vk and extend the function n — > S n to all positive real values by linear interpolation 
between S n and S n+ i. Note that the random function St is Lipschitz with constant (3. 

The following is an immediate consequence of Corollary 5.6. 

Theorem 5.7. Suppose that Vk and S t are as above. If S t and B t are constructed on the 
same probability space (but not necessarily independent), then 

limsup S * - Bt > 1/(4/3). (5.1) 

t^oo logt 

Theorem 2.3.2 of Csorgo and Revesz (1981) says that if the Vk have finite variance 

and 

\S t - B t \ 
hm sup — : = 0, 

t^oo logt 

then the Vk have a standard normal distribution. Our result (5.1) may be interpreted as a 
quantitative version of the same theorem, in the case when \ Vk\ are bounded. A remarkable 
theorem of Komlos, Major and Tusnady (see Csorgo and Revesz (1981) Theorem 2.6.1) 
implies that if the Vk are bounded, then one may construct St and B t on a common 
probability space so that 

limsup-^ — < C < oo. (5.2) 

t^oo log* 

It is striking that one can achieve the same logarithmic order of approximation for a 
Lipschitz function St with independent increments S n — S n -i, as for an arbitrary Lipschitz 
function g(t) with constant f3. Rio (1991) proved that (5.2) holds with C = 9/2 if Vk are 
centered Poisson variables (the estimate had appeared in Section 5 of the preprint; that 
section was not included in the final version of the article, Rio (1994)). No other estimates 
for C seem to be known so (5.1) is our own modest contribution to the field of strong 
approximations. 

6. Open problems. We list a few questions we were not able to answer in this paper. 

(i) Can one prove pathwise uniqueness in Theorem 2.1 if one or both a± and oli belong 
to (-1,0)? 

(ii) Does a result analogous to Theorem 3.7 hold for /?i,/?2 > with /3i — f3 2 < 0? A 
similar question can be asked about the case when (3i < < /3 2 ; in the last case a 
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special solution to (1.1), denned in Lemma 5.2, would have to play an important 
role. Can one generalize Theorem 3.7 to local times corresponding to solutions of 
(1.3) with ct\ and a.2 not necessarily equal to 0? 

(iii) Find the best 7 = 7(0:1, «2, Pi, P2) > in Lemma 4.2. 

(iv) Find the best constants in (5.1) and (5.2). 

(v) Does there exist a unique Lipschitz solution to (2.3) if B t is a fractional Brownian 
motion of index He (1/2,1)? 
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